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ABSTRACT 



Integral field spectroscopy (IFS) is presented for star forming regions in M 33. A 
central area of 300 x 500 pc 2 and the external Hll region IC 132, at a galactocentric 
distance ~ 19' (4.69 kpc) were observed with the Potsdam Multi Aperture Spec- 
trophotometer (PMAS) instrument at the 3.5 m telescope of the Calar Alto Hispano- 
Aleman observatory (CAHA). The spectral coverage goes from 3600 A to 1/im to 
include from [Oll]A3727 A to the near infrared lines required for deriving sulphur 
electron temperature and abundance diagnostics. 

Local conditions within individual Hll regions are presented in the form of emis- 
sion line fluxes and physical conditions for each spatial resolution element (spaxel) 
and for segments with similar Ha surface brightness. 

A clear dichotomy is observed when comparing the central to outer disc Hll 
regions. While the external Hll region has higher electron temperature plus larger 
H/3 equivalent width, size and excitation, the central region has higher extinction 
and metal content. 

The dichotomy extends to the BPT diagnostic diagrams that show two orthog- 
onal broad distributions of points. By comparing with pseudo-3D photoionization 
models we conclude that the bulk observed differences are probably related to a dif- 
ferent ionization parameter and metallicity. 

Wolf-Rayet features are detected in IC 132, and resolved into two concentrations 
whose integrated spectra were used to estimate the characteristic number of WR 
stars. No WR features were detected in the central Hll regions despite their higher 
metallicity. 

Key words: ISM: abundances Hll regions - Galaxies: individual: M 33 



1 INTRODUCTION 

Extragalactic Hll regions provide an excellent laboratory to study star formation processes, evolu- 
tion of massive stars and the properties of the surrounding interstellar medium (ISM). A wealth of 
information can be obtained from the spectral analysis of the bright emission lines and the stellar 
continuum. 

The spectral information from extragalactic Hll regions has been traditionally obtained from 
single aperture or long slit observations of the brightest part of the region. The data is then used 

* jlopez@inaoep.mx 
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to derive the physical conditions of the gas (temperatures and densities) and to estimate abun- 
dances and ionization conditions, as well as characteristics of the ionizing star clusters [masses, 
ages, effective temperatures (T*)]. These result s constitute the main body of our knowledge re 



gardi n g the evolution o f disc galaxies (see e .g. 



1979 



Diaz et al. 



Bresolin et al 



1987 



2004) 



Vflchez et al. 



1988 



Aller 



1942 



Dinerstein 199' 



Searle 



1971 



Skillman 



Smith 



1998 



975; 



Page 



Kennicutt 



et al. 



1998; 



Underpinning this proven methodology to study the emission line spectra lies the tacit as- 
sumption that the observations and the derived measurements are representative of the whole Hn 



region, basically internal variations within the nebula are assum ed to be minima 



tent. The limitation o f this assumption has been long recognized (|Diaz et al 



1987 



or non exis- 



Rubin 



Castaneda et al. 



19891 : 



19921) . however in the majority of the cases any solution has been precluded by 



limitations in the observation technology and optimization of the assigned observing time. Never- 
theless some options have been available to combine spectral information with spatial resolution. 
A simpl e solution is to sweep the Hll region following a direction perpendicular to the slit elon- 



gation (|Kosugi et al 



1995 



Maiz-Apellaniz et alj 120051) in steps separated by one slit width. A 



different approach involves the use of a set of narrow band filters to obtain monochromatic images 
centred on the emission lines that, when combined with neighbouring continuum bands, allows the 
extraction of pure emission images. Tunable filters with Fabry-Perot interferometers offer a more 
flexible scheme by allowing the use of various pass bands, centred at different wavel engths. 



Of particular interest for this work is the technique of integral field spectroscopy (|Allington-Smith 



2006L IFS), also known as 3D or area spectroscopy, in which by the use of integral field unit in- 



struments (IFUs) it is possible to obtain simultaneous spectral and spatial information over the 
observed object. The main advantage of the IFS approach over other techniques is that the ob- 
servation is simultaneous in space and wavelength, producing an homogeneous set of data. The 
drawback of IFS is the small field of view (FOV), thus the selection of the optimal method (tune- 
able filters, long slit sweeping, IFS) for spatially resolved observations depends on the objective 
of each specific project. 

IFS has proved its utili t y in observations of lum i nous and ultraluminqus inf rared galaxies 



(Alonso-Herrero et al. 



2009 



active galaxies (IBarbosa et al. 



2009 



Alonso-Herrero et a 



blue c ompact galaxies (|Cairos et al. 



2009a|b; 



Lagos et al 



2009 



James et al 



2010), the nuc 



Garcia- Marfn et al. 
20091). circ u mnucle ar star forming regions (Dors et al. 



2010 



ei of 



2008) and 



Perez-Montero et al 



201 1|) . Some Hll regions have been studie d with IFS; observations of Hll regions in the outer disc 



of NGC 6946 (Garcia-Benito et al. 



2010L with the same instrument and setup as in this work) 
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among them . IFS o f Hll regions in M 33 has been already obtained (|Monreal-Ibero et al. 



2011 



Relano et al 



2010. for NGC 588 and NGC 595 respectively), in both cases the PMAS (Pots- 



dam MultiAperture Spectrophotometer) instrument at the 3.5m telescope of the Observatorio As- 
tronomico Hispano-Aleman in Calar Alto (CAHA) was used. They compared the distribution maps 
for different parameters against the reported values in the literature, obtaining a good agreement 
when the IFS data is integrated in a single spectrum and noting that for some quantities (e.g. ex- 
tinction) to assume a single value for the whole region is not strictly correct. For some of the 
empirical abundance estimators a spatial dispersion is present, however the variation is within the 
estimated observational error and in such case a uniform distribution may be assumed. 

The general objectives of this work are, using IFS data for the centre and an external region of 
M 33, to characterize the extent of the internal variations in Hll regions for measured and derived 
values; to test whether the description of the whole Hll region with a single value for different 
parameters is valid; to map the massive star content using characteristic spectral features (i. e. 
Wolf Rayet stars) and to compare internal variations in high metallicity (inner) vs. low metallicity 
(outer) Hll regions. Also, given that the central region remains an ill defined zone with scarce 
spectroscopic observations, IFS may contribute to its better characterization. 

The organization of the paper is as follows. The observations, data reduction and integrated 
emission line maps are described in Section 2. Section 3 describes the physical conditions of 
the gas. Section 4 concerns the chemical abundance determinations. Section 5 is devoted to the 
discussion of diagnostic diagrams. The detection of Wolf-Rayet features is presented in section 6. 
Analysis in projected shells is shown in Section 7. Finally the conclusions are presented in Section 
8. 



2 OBSERVATIONS AND DATA REDUCTION 



2.1 Object selection 



M 33 is the third brightest member of the Local Group (|Van den Bergh Il2000|) . Given its prox- 
imity, large angular size and low inclination, it is an ideal candidate for IFS observations, having 
a rich content of Hll regions in the central region and across the disc at different galactocentric 
distances allowing the exploration of recent star formation a ctivity in a wide r ange of physical 



conditions. Hundreds of nebulae have been charted in M 33 (IHodee et al 



abundance and gradient studies (|Vflchez et al. 



1988 



Crockett et al 



2006 



2002) and included in 



Rosolowskv & Simon 
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Table 1. M 33 properties. 



Designations Messier 033, Triangulum Galaxy, 

NGC 0598,UGC 01117, PGC 0058 U 
Classification SA(s)cd Hll 



Major Diameter 70.8 Arcmin 

Minor diameter 41.7 arcmin 

Position Angle 23 deg 

Distance 840 kpc (Freedman et al. 2001) 

Redshift -0.000597 ± 0.000010 

Inclination angle 53.52 deg (Corbelli & Salucci 2000) 

PPak scale 4.07 pc/arcsec (10.91 pc/fibre ) 

Source is NED unless otherwise specified. 



2008 



Bresolin et al.ll2010h . Tabled] summarizes the main properties of the galaxy. Values are ob- 



tained from NED0 wherever no references are given. 

As a first step of a larger project aimed to trace the chemical abundances, gradient and dis- 
persion of physical conditions with high spatial resolution across M 33, we obtained IFS data 
of the centre and the external Hll region IC 132, located to the NW of the centre of M 33 at a 
galactocentric distance of ~ 19' or 4.69 kpc (see Fig. [I]). 

The observed locations represent two scenarios of star formation in spiral galaxies with ex- 
treme different environments: outer spiral disc and central regions. Notorious differences between 
the two (disc vs central) exist regarding the star formation rate (20 vs 1000 M Q yea r" 1 ), star forma 



tion t i mescales ( lMyr vs 1 Gyr), electron density (10 vs 10 cm 6 ) among others (jKennicutt et al. 



1989; 



Kennicutt 



1998). 



2.1.1 The central region ofM 33 

The cent ral region o f M 33 has been the subject of various studies given its notorious star formation 
activity (|Keelll 19830 . From the fact that the dominant stella r population in the ce ntre of M 33 is 
young (< 1 Gyr) and metal rich (up to solar metallicities), IVan den Berghl (| 199 II) suggested that 
the stars are the product of intense star forming periods originated by inflow of gas. 

Star formation in the centre of spiral galaxies has been associated with the existence of bulges 
and bars. However the presence of a bulge in M 33 is controversial. The excess in surface bright- 
ness in the central region of M 33 (|Kentlll987|) is n ot necessarily due to the existence of a bulge 



but may be explained by a large scale diffuse halo (Bothun 



By using H-band photometry of individual stars 



Minniti et al. 



1992, based on 12 jim photometry). 



19931) concluded that a bulge exists 



and that it has experienced a period of star formation less than 1 Gyr ago. The existence of a bar 

1 This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA). 
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in M 33 is also an open issue, iRegan & Vogell (I1994|) propose that if a bar exists it has a length 
of 3' with a position angle of 90°. They propose also that the bar-like emission could be due to 
the continuation of the spiral arms. However they claim that there is a lack of evidence in order 
to decide in favour of one or the other explanation. In any case, in M 33 the bar-like structure 
contributes only 4% to the total light out to 1'.5, in contrast with a 40% contribution in strong 



barred gal axies (jB lackman 



19831). Also the existence of a massive central black hole has been 



discarded (K ormendv & McClure 



1993 



Massev et al. 



1996) and no AGN-like emission line ratios 



have been reported. However the centre of M 33 harbours M 33 X-8, the brightest ultraluminous 
X-ray source (ULX) in the local group. The observational evidence indicates that M 33 X-8 emis- 
sion originates in a stellar-mass black hole (M Rff ~ 1OM0) w ith super- Eddington accretion rate 



(IFoschini et al. 



2004; 



Weng et al 



2009; 



Middleton et al. 



20111) . M 33 X-8 position which coin- 



cides with the photometric centre of the galaxy, is indicated with a cross symbol in the left panel 

of Fig. m 

The central region is an excellent scenario to study the star formation process in relatively 
high metallicity and density and low excitation environments, offering a challenging scenario de- 
viating from the conditions of the majority of Hn region studies. Such conditions represent a 
test to the different abundance calibrators, generally based on low metallicity Hn region samples 



(IKinkel & Rosa 



1994; 



Bresolin et al 



20041) and extrapolated to the high metallicity regime. 



2.1.2 The outer disc region IC 1 32 

IC 132 is one of the most external Hn regions in M 33, located 4.69 kpc (~ 19') away from the 



centre. It ranks amon g the 



and g radient studies (Allei 



1977 



brightest in the galaxy and 



1942 



Searle 



1971 



Smith 



las been included in chemical composition 



1975 



K witter & Aller 



1981 



Magrini et al. 



2007 a|). Wolf-Rayet features we re detected from narrow band interference filter surveys (IBoksenberg et al 



D'Odorico & Benvenuti 



19831) . Metallicity has been reported ranging fr om 12+log(0/H)=7.85 



to 12+log(O/H)=8.08 with typical errors ~ 0.05 (Ma grini et al.ll2007b . 



2010). 



IC 132 shows the high excitation, low density and low abundance environment typical of disc 
regions, for which most of the Hn studies, physical conditions and abundance estimators have 
been developed. 
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Figure 1. Position of the observed fields over M 33. Hexagons are 74" in diameter. The diamond marks the centre of the galaxy with coordinates 
RA=01 h 33 m 50".89, DEC=+30°39'36."80 (J2000). The outer region IC 132 has coordinates RA=01 h 33 m 15 s .90, DEC=+30°56'44".00 (J2000). 
North is up and east is left. 



2.2 Observations 



Observations were obt ained with the CA HA 3.5m telescope, using the PMAS instrument in fibre 



Package (PPak) mode (IRoth et al. 



2005). With this configuration the IFU features 331 fibres for 



the science object packed in an hexagonal array with a FOV of 74" x64", surrounded by 6 bundles 
of 6 fibres each at a distance of 72" from the centre to sample the sky. 15 extra fibres, not on 
the focal plane, can be illuminated separately for calibration purposes. Although the unit contains 
more fibres to provide protection against mechanical stress, the total number of active fibres adds 
to 382. Each fibre has a projected diameter of 2". 68 in the sky or 10.91 pc at the assumed distance 
of M 33. The filling factor of the science packet is 60%. The detector use d was a thinned, blue 



enhanced, 2Kx4K CCD (SITe ST002A) with 15/xm pixels (IKelz et al. 



2006). 



The centre of the galaxy and the external region IC 132 w ere observed between September 
7th and 9th 2007. The position of the fields over an Ha image ( Cheng et al. 1 1997 ) are in Fig. [T] 
Detailed identification is in Fig.El the instrument FOV is overplotted on NOAO Science Archive 



images (|Massey et al. 



2006) taken with a 50 A FWHM filter centred on Ha. The fields for the 



centre were overlaid by one fibre width. 
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Figure 2. Identification of the observed fields. Images are narrow filters centred in Ho from the NOAO Science Archive l lMassev et alj 2006). 
Overlaid is the PMAS-PPak hexagonal FOV. For the centre of M 33 (left) two consecutive fields are overlapped by one fibre line. The cross in the 
central regions indicates the centre of the galaxy coincident with the ultraluminous X-ray source M 33 X-8; for IC 132 (right) the cross marks the 
centre of the region. The scale in parsecs corresponds to the assumed distance to M 33 of 840 kpc. 

Observations were performed at two different grating angles (GROT) with the 300 grooves 
mm (V300) grating. The value GROT=-75 was selected to cover from 3591 to 6996 A (op- 
tical) with a dispersion of 3.4 A pixel -1 , and GROT=-72 to cover from 6873 to 10186 A (near 
infrared) and a dispersion of 3.2 A pixel -1 . In addition high spectral resolution spectra were ob- 
tained using the V1200 grating at GROT=-55.5, covering from 6100 to 6650 A with a dispersion 
of 0.64 A pixel -1 . In all cases the CCD was sampled in a 2x2 binning mode and the values quoted 
above are the effective ones after the resampling. A journal of observations is shown in table [2] 
The stars BD+17 4708 and BD+28 4211 were observed for flux calibration in the optical and 
NIR ranges respectively. Various continuum illuminated and HeHgCs+ThAr lamp frames were 
acquired, interleaved and as close as possible to the science observations to account for flexures in 
the telescope. Sky flats were acquired at the beginning of the night and biases at the end. For each 
object 3 dithered pointings were taken in order to fully cover the instrument FOV. The positions 
are D1=(0,0), D2=(l. 56,0.78) and D3=(l. 56,-0.78), where the displacements are in arc seconds 
relative to the initial position Dl. The average seeing was 1.0"; higher values as well as the pres- 
ence of clouds are indicated in the last column of table |2] Half of the third night was lost due to 
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Table 2. Journal of observations. 



Object 


Spectral range 


Dispersion 


Exposure Time 


Median 






[A] 


[A pix- 1 ] 


[s] 


Airmass 


Notes 


centre PI Dl 


3591-6996 


3.4 


500 x 3 




seeing ~1".46 


centre PI D2 


3591-6996 


3.4 


500 x 3 


1.49 


seeing ~1".43 


centre PI D3 


3591-6996 


3.4 


500 x 3 






centre P2 Dl 


3591-6996 


3.4 


500 X 3 






centre P2 D2 


3591-6996 


3.4 


500 x 3 


1.01 




centre P2 D3 


3591-6996 


3.4 


500 x 3 






IC 132 Dl 


3591-6996 


3.4 


500 X 3 






IC 132 D2 


3591-6996 


3.4 


500 X 3 


1.02 




IC 132 D3 


3591-6996 


3.4 


500 X 3 






centre PI Dl 


6873-10186 


3.2 


500 + 700 x 2 




seeing ~1".33 


centre PI D2 


6873-10186 


3.2 


700 x 3 


1.49 




centre PI D3 


6873-10186 


3.2 


700 X 3 






centre P2 Dl 


6873-10186 


3.2 


700 x 3 






centre P2 D2 


6873-10186 


3.2 


700 X 3 


1.04 


clouded, seeing ~ 1 " .48 


centre P2 D3 


6873-10186 


3.2 


700 X 3 






IC 132 Dl 


6873-10186 


3.2 


500 X 3 




clouded 


IC 132 D2 


6873-10186 


3.2 


500 X 3 


1.07 


clouded 


IC 132 D3 


6873-10186 


3.2 


500 X 3 




clouded 


centre PI Dl 


6100-6650 


0.64 


500 X 3 






centre PI D2 


6100-6650 


0.64 


500 X 3 


1.19 




centre PI D3 


6100-6650 


0.64 


500 X 3 




guiding stopped momentarily 


centre P2 Dl 


6100-6650 


0.64 


500 X 3 






centre P2 D2 


6100-6650 


0.64 


500 X 3 


1.01 




centre P2 D3 


6100-6650 


0.64 


500 X 3 






IC 132 Dl 


6100-6650 


0.64 


500 X 3 




clouded, seeing 1".50 


IC 132 D2 


6100-6650 


0.64 


500 X 3 


1.04 


clouded, seeing 1".40 


IC 132 D3 


6100-6650 


0.64 


500 




clouded, seeing 1".60 



The median seeing was 1".0, those pointings with seeing above 1".3 are indicated as well as other non photometric conditions. 
The median airmass is given for each 3 dithering positions. For the third night no bias or sky flat frames were taken. 



weather which limited severely the observing time and neither bias nor sky frames were taken. In 
this case the bias was determined from the overscan region of each frame. 
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Figure 3. Integrated visual and near-IR spectra for IC 132. Relevant emission lines are identified. See text regarding special considerations for 
[OlII]A4363A and Garci'a-Benito et al. 2010 for the treatment of the near IR [Sill] lines. 



2.3 Data reduction 



The data reduction is based o n the set of routines and tools group ed in the E3D (|Sanchez 



Sanchez. Becker & Kelz 



20041) and R3D ([Sanchez & Cardiel 



2004; 



2005) packages. These were com- 
plemented with visualization and analysis tools developed by us, together with standard IRAFc 
utilities. 

The pipeline includes spectra tracing, extraction, fiber flat correction, flux calibration, sky sub- 
traction and datacube building. Once the 3D datacube was assembled, it was additionally corrected 
for differential atmospheric refraction, atmospheric correction for the fSllll ne ar infrared lines and 



Garcfa-Benito et al 



(|2010f) and in the Comple- 



extinction. Details of the reduction can be found in i 
mentary Online Material. 

The final integrated spectra for the IC 132 region is in figure |3] with relevant lines identified in 
the optical and NIR spectral ranges. 



2.4 Emission Line Maps 

Emission lines were fitted with a gaussian, and continuum bands at both sides of the emission were 
measured additionally to estimate the error in the continuum placement. Given the low resolution 
of the data (FWHM effective resolution Rfwhm - 700 at A = 5500 A and Rfwhm - 1 100 at A 

2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in 
Astronomy (AURA) under cooperative agreement with the National Science Foundation. 
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= 8500 A) for the first and second nights, some of the lines had to be deblended by fitting multiple 
gaussians. This was the case for the triple blend of Ha and the two [Nil] lines at AA 6548,6584 A 
and the [Sll] AA 6717,6731 A doublet. In each case a single gaussian is assigned to each line and 
the best solution is determined using a non-linear least-squares fit. For the third night, the higher 
resolution (Rfwhm — 3100 at A = 6500 A) allows the separation of the Ha-[Nll]AA 6548,6584 A 
system as well as the measurement of the faint [Sill] A6312 line, at least for the brightest spaxels 
ofIC 132. 

Examples of 2D maps are shown in figure |4]for IC 132: 1(H/3), H/3 equivalent width [EW(H/3)] 
and I([Om]A 4363), and in figure |5]for the central zone: I(H/3), EW(H/3) and I([Sm]9069). The 
line maps have S/N > 2.0. 
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3 PHYSICAL CONDITIONS OF THE GAS 

Physical parameters describing the ionized gas were derived using the line fluxes measured in 
12.41 To obtain these parameters we used the TEMDEN routine included in the STSDAS nebular 



package for IRAF. This task so 



imation (Shaw & Dufoui 



ves the statistical equilibrium equation for the 5 -level atom approx- 



19951). The transition probabilities have been updated with recent v alues 



computed for the [On], [Om], [Nil], [Sii] and [Sin] io ns by Fisch er & Tachiev 



collisi on strength values were taken form 



[OIII 



Hudson & Bell 



(120051 ) for [Nil], 



Taval 



d2007h for TO 



I], 



(2004 



Aggarwal & Keenanl ( 



Ramsbottom & Belli < 19961) for [Sii] and 



2006). The 



1999) for 



Tayal & Gupta 



(119991) for [Sm]. 

Computations were made for each individual spaxel where the required lines are detected with 
enough S/N and also for the shells analysis described in section SJ7J The S/N > 5 requirement 
was met by the intense lines ([Olll]A5007, Ha, H/3) for all the spaxels, while for the auroral lines 
[Om]A4363 and [Sm]A6312, S/N > 2 was reached only in the central part of IC 132. 

The electron density (n e ) was estimated from the [Sll]A6717/[Sll]A6731 line ratio, and for the 
electron temperature (T e ) we used the faint auroral [Om]A4363 and [Sill] A6312 lines, allowing the 
use of the [Oni]AA4959,5007/[Oiii]A4363 and [Siii]AA9069,9532/[Sm]A6312 ratios to compute 
T e ([Olll]) and T c (Slll) respectively. Due o strong atmospheric absorptions around 9300 A, we 
defined the [SIII] 9532 line by its theoretical ratio to the 9069 A line. 



3.1 Electron densities 

To esti mate n P we used the ratio RS2 = [SII]A6717/[SII]A 6731 and the TEMDEN routine in 



IRAF (Shaw & Dufour 



19951) . Fig. [6] shows the RS2 ratio for the external region IC 132 and for 
the centre of M33. For IC 132 the RS2 distribution is quite uniform with a mean value around 
1.45 and although the RS2 ratio shows variations in the brightest parts, the conversion to n c would 
translate into variations within the low density regime, which is of no significant impact in further 
parameter estimations. 

The central region of M33 shows a broad distribution with the brightest Ha region having RS2 
values around 1.4 indicating n c values also in the low density regime, while the northern much 
fainter Ha extension shows some spaxels with unrealistically large RS2 ratios. From the RS2 ratio 
it seems that a structure exists and in such case the hypothesis of a uniform value across the nebula 
would be too simplistic. 

From the n c maps it is obvious that many pixels show an RS2 ratio above the theoretical 
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value for low density limit and that it decreases down to the lower ratio limit, indicating that the 
region spans the whole range of n c from very high to very low densities. There is positively some 
uncertainty from deblending the [Sn] lines, and this is not exclusive of our data. Previous RS2 
determinations from IFS data also show ratios wit h significant excursions above the theoretical 



limit, from data obtained with the same instrument (IRelano et al.ll2010; 



vlonreal-Ibero et al. 



2011) 



2009), even with the [Sn] 



or other integral field units (IFU) in a different telescope (|Lagos et al. 
lines resolved. 

However the problem does not seem to be exclusive of IFS data. From a literature compilation 
and their own observations, all obtained with long slit spectroscopy, Kennicutt et al. (1989, figure 
5) reports an RS2 ratio peaked at 1 .3 and with a significant number of Hn regions, mainly in spiral 
disks, above the 1.4 theoretical limit. Also from long slit spectra, Zaritsky et al. (1994, figure 4) 
report a similar behaviour for RS2. 

When RS2 ^ 1 .4 is obtained, it may be assumed that n c is below 10 cm~ 3 however, the assign- 
ment of a true density is very uncertain (it begins to be uncertain as RS2 reaches the asymptotic 
limit even before 1.4). Although the safe way to proceed in such cases is to assume n c =100 cm" 3 

- in any case for temperature and abundance determinations the density plays a second order role 

- it seems that the theoretical ground for the determination of n e may need an adjustment to the 
atomic data and further investigation in the topic is required which is beyond the scope of this 
work. 

This problem, that seems to be affecting both ends of the theoretical calibration, could be traced 
to uncertainties in the determination of the atomic data used, specifically the collision strengths for 



the low-end and the transition proba bilities for the high-end, not only for fSlll but for 



cato rs based on different line ratios (Mc Laughlin & B ell 



1998 



density indi- 



Copetti& Writzlll2002 : 



Wang et al 



20041 Gary Ferland, private communication). We concur with Wang and collaborators, that a full 



independent study with ad-hoc data to validate previous results will be crucial. 



3.2 Electron temperatures 

We have determined the T c ([Om]) and T c ([Sm]) electron temperatures using the direct method 
based on the line ratios [Om]AA4959,5007/A4363 and [Sm]AA9069,9532/A6312 respectively. 

Only in the brightest parts of the core of IC 132 it was possible to derive reliable T c ([Om]) 
and T c ([Sm]). Unfortunately the sky at Calar Alto suffers from light pollution produced by street 
lamps, and among the several emissions in the visible range is the Hg 1 4358.3 line, which is strong 
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Figure 6. [SII]A6717/[SII]A6731 maps for IC 132 (left) and for the central zone of M33 (right). Notice the different vertical scale in the histograms. 



enough (iSanchez et al.ll2007l) to limit the use of the [Om]A4363.2 line in the temperature determi- 
nation by the direct method in low redshift systems. This is particularly bad for M33 because of 
its blue-shifted spectrum (z=-0. 000597) that brings the [Olll] line to A=4360.61, close enough for 
it to blend with the Hgl line. 

Nevertheless, by performing a careful sky subtraction, we were able to obtain upper limits 
for the [Om]4363 line. This upper limit was used to compute an upper limit of T c ([Om]) and to 
estimate the Ionization Correction Factor (ICF) for S. T c ([Olll]) was computed using the ratio R o3 
= ([Oiii]4959+[Oiii]5007)/[Oiii]4363 and the result is shown in figure [7] top left panel. 

For T e ([Oii]) we face the problem that while the [Oil] A A 3727,3729A doublet was detected, 
the [Oll]A7320 and A7330 lines were not. Thus we have to relay on results from fits to photoion- 



ization models an d their comparison with data with direct T e ([On]) and T e ([Om 



Diaz et al. 



) determinations 



Stasihska ( 



1990) 



(I2007h assume T([On]) ~ T([Sm]), while Izotov et al. 2006, using 
models, suggest taking T c ([0n]) = T c ([Olll]) if no other T e ([0n]) determination is possible, how- 
ever it is pointed out that in such case, the error based in photoionization models with Z=Z Q 
is equivalent to underestimating [Oil] AA3727,29 by about 40-50%. In some cases T c ([Oil]) = 
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T e ([Nll]), however to estimate T ([Nn]), the weak [Nil] A5755 line is required and in this case 
also an underestimation of [Oil] is obtained (IHagele et al] |2008). In any case, [Nil] A5755 is not 
detected in our observations. 

It seems that no strong correlation exists between the two temperatures in Hn regions or in Hn 
galaxies and that the effect of meta llicity, density and maybe even other factors should be taken 



into account. 



Kennicutt et al. 



l2003i ) find only "a hint of correlati on", using their observations and 



data from the literature. N evertheless, they use the lGarnettl (11992) expression to estimate T([Oll]). 



Hagele et al. 



2006, 



2008|) also show that the cor relation between T P ([Oii 



weak and favour the use of the density dependent 
Here we apply the widely used relation from 



Perez-Montero & Diaz (2003) relation. 



Garnett ( 



) and T e ([Om]) is very 



1992) and the more recent fits from 



Perez-Montero & Diazl (|2003|) to estimate T e ([Oil]) from the upper limit of T e ([Olll]). Fig. |7] 



shows the results from both estimates in the top middle and right panels. For T c ([Oii])pmo3 
n e =100 was used for the whole region. 



The electron temperat ure determined using th e 



Garnett 



about 1600 K higher than 
that for n„=100 an d T P (TO llll) > 15000, the 



below the 



(119921) method gives a mean T ([Oil]) 



Perez-Montero & Dfaa (120031). This difference may be traced to the fact 



Perez-Montero & Diaa (|2003|) model is systematically 



Garnettl(l 19921) relation. If n c =10 is assumed, then the T c ([0n]) would increase at least 



2000K, bringing the two methods to better agreement. 
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Figure 7. Top left to right: r e ([O III]) up per l imit for IC 132. T e ([OlI]) estimated from the relations provided by iGarnetdd 1992b andjPerez-Montero & Diaz (2003). Bottom left to right: T c ([SlII]) obtained from the 
direct method, using Hagele et al. (200eJ) and lDi'az et all fc007l) . Here and in following figures, the histogram of the distribution also indicates the total number of pixels used (N), the mean (m and solid blue line) 



and the standard deviation (s and red dashed lines). 
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The sulphur electron temperature T c ([Slll]) was computed for spaxels where [Sm]A6312 was 
detected with S/N > 2, using the ratio R 5 23=([Siii]9069+9532)/[Sm]6312. The intensity of [Slll]9532 
was obtained from the theoretical relation [Sm]9532 = 2.44 x [Sm]9069, as the A9532 A line is 
more affected in these spectra by telluric emissions and a bsorptions. 

T e ([Slll]) was also estimated using the fit obtained by (IHagele et al.ll2006l H06) from a large set 
of observations that include Hn galaxies, Giant Extragalactic Hn regions and diffuse Hn regions 
in the Galaxy and the Magellanic clouds, given as T e ([Sm]) = 1.19T e ([Olll]) - 0.32 

The set includes objects with T c ([Sm]) up to 24,000 K. This extends the relation to higher 
temperatures but also increases the uncertainty in the fit given that at the high temperature regime 
fewer objects exist and they exhibit large errors. In this case we remark that T c ([Om]) is an upper 
limit and caution should be taken for i ts interpretation. 



We also used the SO23 parameter (|Dfaz & Perez-Montero 



2000; 



Diaz et al. 



2007), defined as 



S 23 _ [Sn]AA6717, 6731+ [Slll]AA9069, 9532 
23 ~ Chz ~ [On]AA3727, 3729 + [Om] AA4959, 5007 ( } 

which has weak dependence on the ionization parameter. T c is then obtained as T e ([Slll]) = 0.596 

- 0.283(log S0 23 ) + 0.199(log S0 23 ) 2 . 

Although this calibration was obtained for high metallicity objects and its use is justified for 
the central region, we also applied it for IC 132 given that the SO23 values obtained are within the 
parameter validity range of the original fit. 

Fig. [7]maps the distribution of T c ([Slll]) computed using H06 and SO23 methods for IC 132. 
In the first case the dependency on T e ([Om]) and consequently on [Om]A 4363 is reflected in the 
facts that only the area where [Om] is detected with S/N > 2 has T e ([Sm]) and that these values 
are higher than the ones obtained with the other methods, being T c ([Olll]) an upper limit. 

Besides the fact that SO23 has low sensitivity to the ionization parameter, the difference in 
scatter may also be traced to the difference in errors associated with the lines used. For the H06 
method the large error of the [Olll]A 4363 line is propagated to give a standard deviation of about 
3800 K in T c ([Sm]), while for S0 23 all the involved lines having S/N > 10, give a standard 
deviation of just 700 K in T c ([Sm]). The low mean T e ([Slll])=9800 may have its root in the fact 
that the SO23 calibration was obtained mainly with high metallicity objects. 

Fig. [8]shows the temperature distribution for the central region of M33 using the SO23 method. 
It gives also a very uniform temperature for all the measured spaxels in the field. 

If the scatter in temperature can be mainly associated to the ionization structure and the SO23 
effectively cancels it then it would be worth to investigate the inclusion of more low metallicity 
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Figure 8. T c ([SlII]) obtained as in lDfaz et all J2007l) for the centre of M33. 



objects in the T c ([Sm]) SO23 fit. This m ay modify the re lation to produce higher temperatures in 



the low metallicity range [see figure 9 of 



Diaz et al. 



2007)] 



4 CHEMICAL ABUNDANCES 



4.1 Direct method 
Oxygen 

Taking the upper limit temperatures obtained for the central zone of IC 132, we computed the 
lower limit ionic abundances using the IONIC task from the N EBULAR package for I RAF, also 



based on the 5-level atom statistical equilibrium approximation (|De Robertis et al. 



1987). This task 



requires additionally the electron density n c and the emission line flux relative to H/3 of an atom 
in specific ionization stages. For ++ /H + we used the uppe r limit T P (rOllll) while for Q + /H + 



we adopted the T c ([0n]) obtained from the relations given by 
total oxygen abundance was obtained as 
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O 
H 



o- 



+ 



Or 



(2) 



H+ H+ 

and is displayed, together with the ionic values, at the top of figure [9] 
Sulphur 

Sulphur ionic abundances were obtained using the expressions derived by 
fits to results from the IONIC task. The unobserved S 3+ can represent an important contribution in 
high ionization zones, thus in order to estimate the total sulphur abund ance an ionization correction 



Diaz et al. 



(2007) from 



Barker! (11980h . Although 0+ and C++ are 



factor (ICF) must be included. We used the ICF from 
lower limits, their inclusion in the sulphur ICF yields a better S estimation than just ignoring the 
correction factor. 

The bottom of figure|9]shows the ionic and total sulphur abundances for IC 132. For the central 
region the corresponding maps are in figure [TO] The spatial distibution in IC 132 seems to have a 
wider dynamical range than the central zone in S + and S ++ . This variation is primarily produced 
by a diagonal zone in the lower left corner of IC 132, which we call a "wall" that could probably 
be attributed to a zone of diffuse radiation (see figure ??). 
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4.2 Empirical methods 



The direct method is useful in a restricted range of abundances. For systems with oxygen abun- 
dance larger than ~ O.2(O/H) , as the metallicity increases, so does the cooling by IR fine structure 
lines and the rOlHlA4363A line becomes extremely weak, requiring high S/N observations to be 
detected. IIKennicutt et alj ( 120031) estimate oxygen abundances ~ O.5(O/H) as the practical upper 
limit to detect [Om]A4363A]. 

The difficulty of direct abundance estimation has led to the development of empirical or strong 



line methods amply used in the literature. 



4.2.1 The i?23 parameter 



The empirica l parameter R 93 =( rOiil3727,29+rOml4959, 5007)/H/3, also known as 23 , was first 
proposed by IPagel et al.l (|1979|) and is widely used as an abundance indicator. Several calibra- 
tions exist for i? 23 , revised in the light of new observations with direct T c determinations or im- 
provements to photoionization models. Some reformulations of the parameter have also been pro- 
duced in an effort to compensate for the effec t of the stellar temperature and ionization parameter 



(Edmunds & Pagel 



1984; 



Zaritskv et al. 



1994 



Kobulnickv et al. 



1999 



Pilyugin 



200 lib) 



R23 has two main drawbacks. Firstly is the fact that the metallicity mapping is double valued. 
For each value of R 2 3 there is a low metallicity value corresponding to the lower branch of the 
relation and a high metallicity estimate corresponding to the upper branch. Besid es, the transition 



Kennicutt et al. 



betw een branches depends on the calibration and is ill defined [see e.g. Fig. 12 of 
(|2003|) 1. Secondly, most of the data points tend to lay close to the knee of the relation, i.e. in the 
region in which there is no dependence of the R 2 3 parameter on metallicity. The R 23 method is in 
practice a reliable abundance estimator only for the metal rich or very metal poor systems. 

Here we used the calibratio n from McGaugh ( 199l|) . based on Hn region photoionization 



1998) and taking into account the effect of the ion- 



models obtained with Cloudy (IFerland et al. 
ization parameter . The analytic expressions for the upper and lower branches are taken from 



Kobulnickv et al. 



dl999t) . The estimated uncertainty for this calibration ranges from 0.1 to 0.2 
dex when the lines are detected with S/N > 8, reaching up to ±0.25 dex in the turnover zone 
around 12 +log(0/H) ~ 8.4. 

To break the R 23 degeneracy the [Nll]A6584/[0n]A3727 ratio is used as an initial guess of 
metallicity helping in the selection of the appropriate R 23 branch (McCall, Rybski & Shields 
1985). This ratio shows a weak dependence on the ionization parameter and a strong correla- 
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tion with metallicity (|Kewley & Do pita 



2002). Anot 



proper R23 bra n ch is [Nll]A6584/Ha: (IRaimann et al. 



ler di agnostic ratio used to distinguish the 



Denicolo et al 



2000) the N2 index, calibrated in [O/H] by 



(12002|) . However this ratio has a strong dependence on the ionization parameter 



(see below). 



Fig. dU shows R 23 as a function of the [Nn]A6584/[0ll]A3727 ratio. The Kewley & Ellison 



(|2008|) empirical separation between the upper and lower branches of R23 at [Nn]A6584/[0ll]AA 



3727 = -1.2 [12 + log (O/H) ~ 8.4] is used. For IC 132 the majority of the pixels fall within the 
lower branch regime while for the central zone the opposite is observed. The central region seems 
to follow a continuous distribution of R23, while in IC 132 most of the points are clustered between 
log(R 2 3) ~ -1.2 and ~ 0.8 an d the distribution shows a particularly abrupt cutoff at log(R 2 3) ~ 0.8. 

The empirical division by llCewley & Dopita ( 2002 ) was obtained combining observations and 
photoionization model grids. These models show a dispersion below [Nll]A 658 4/[OlllA 3727 

.2 f or different values of ionization parameters and metallicities [Fig. 8 from Kewley & Dopita 
(I2002|) 1. However the correlation is remarkably extended up to log(R 2 3) > 1.6 for the ce ntre. This 



tail is produced by the faint (diffuse) radiation, not included in lKewley & Ellison! (120081) models. 

The resulting maps for the R23 estimator are shown in the left panels of Fig. [T2]for IC 132 and 
figure[[3]for the centre of M33. 

Another empirical indicator based on the R23 method is proposed by IPilyuginl (|20Q 1 allbllch 
where the excitation parameter P is introduced to compensate for R 2 3 variations along the region 
produced by differences in the ionization parameter. The upper branch applies to 12+log(0/H) 
> 8.2 and the lower branch to 12+log(0/H) < 8.2 Results from an estimation spaxel by spaxel 
are shown in Figs. [T2]and [13] for IC 132 and the centre respectively. In this case we used the 



[Nn]A6584/[0ll]A3727 ratio to select the branch and break the R23 degeneracy. 



4.2.2 The N2 parameter 
The N2 index, defined b y 



Denicolo et al. 



(120021) as N2=[Nll]6584/Ha, was previously used by 



Storchi-Bergmann et al.l (|1994|) as an empirical abundance estimator for star forming galaxies. Its 
usefulness is based on the fact that it is easy to observe up to relatively high redshifts and, due to 
the closeness of the lines involved, it is unaffected by red dening and can even be measured directly 
from uncalibrated spectra. Here we use the N2 fit from [Denicolo et al.l (|2002r) . given as 12 + log 
(O/H)=9.12+0.73N2. The resulting maps from the N2 estimator are shown in the lower panels of 
Fig. OH for IC 132 and in figured!] for the centre of M33. 
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Figure 11. [NII]/[OII] vs R23 to select the upper or lower branch for each pixel. Left: IC 132; right: central zone. The dividing dashed line is from 
iKewlev & Ellisonl {2OO8) calibration. Only those pixels with errors lower than 0. 1 dex in each axis are plotted. 

It is remarkable that the abundance based on N2, shows practically no variation in the maps. 
One may conclude from this indicator that the oxygen abundance has one single value for the 
whole region. It seems that all the variation within the region can be traced to the error which 
effectively increases with radius b oth for IC 1 32 and for th e bright central condensation hereinafter 



called BCLMP 93 as identified in 



Boulesteix et al. 



1974). 



4.2.3 The 03N2 parameter 



19791) as an abundance estimator, 



Pettini & Pagell (I2004h calibrated the 03N2 ratio (lAUoin et al. 
giving the relation 12 + log (0/H)=8.73 - 0.32 x 03N2 where 03N2 = log([Om] A5007/H/3)/([Nll] 
A6584/Ha). 

The resulting maps of 12 + log(0/H) from the 03N2 estimator are shown in the right panels 
of figure d2]for IC 132 and figure [T3]for the centre of M33. 

The distribution of O/H obtained via N2 and 03N2 are very similar; the scatter is smaller when 
the abundance is derived using 03N2 than when using N2 in IC 132. 



4.2.4 The S23 parameter 

The S23 parameter, proposed by 



Vflchez & Esteban 



(11996 ) that uses the strong [Sn] and [Sin] 



lin es, in analogy to R , ? is de fined as S 23 =([Sn]6717, 6731 + [Sm]9069, 9532)/H/3. 



Christensen et al. 



(|1997|) calibrated S23 as a sulphur abundance indicato r based on observa- 
tional data compiled from the literature, while Diaz & Perez-Montero ( 200ol) proposed S23 as an 
alternative to R 23 as oxygen abundance indicator. S23 presents some advantages over R23: it is 
single valued up to solar metallicities and has less dependence on the ionization parameter and on 
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Figure 12. Empirical metallicity estimators for IC 132. Top: R23 with the calibrations from McGaugh ( 1991) and Pilyugin (2001c). Bottom: N2 
and 03N2. 

the effective temperature of the ionizing stars. On the other hand, it needs wider spectral coverage 
to detect the near-IR [Sill] lines and corrections for the unseen S 3+ i on. 
We estimated the sulphur abundance (S/ H) using the relation from 



Perez-Montero et al 



Pere z-Montero & Di az (2005). Both are fits obtained 



(2006), 



while for (O/H) we used the expression from 
from observed data sets. 

Sulphur and oxygen abundance maps obtained from the S23 indicator are shown in figures [14] 
and Q3] respectively. The mean 12+log(S/H) 52 3 value for the central region is 7.02 with a standard 
deviation of 0.34 while for IC 132 it is 9.5 with a smaller standard deviation of 0.28. The estimates 
of S/H and O/H for IC132 using S23 show little scatter, albeit slightly larger than that for the O/H 
abundance using the N2 or 03N2 estimators. 
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Figure 13. Empirical metallicity estimators for the centre of M33. Top: 
Bottom: N2 and 03N2. 
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R23 with the calibrations from McGaugh jl99ll) and IPirvu girl bOOld) . 
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Figure 14. 12 + log(S/H) abundance using the S23 indicator for IC 132 (left) and for the centre of M33 (right). 



4.2.5 Ar 3 3 and S3O3 calibrations 



Stasiriska 



Based on a compilation of Hn regions in spiral galaxies and blue compact galaxies, 
J2OO6 ) proposed the use of two line ratios (Ar 3 3 = [Arm]A7135/[Oii]A5007 and S 3 3 = 
rSml A9069/rOll ll A5007) as metallicity indicators. The calibration of the relation is based on the 
Pilyuginl ( 2001a ) method. In $7] we will analyze the results separating the data in sections of equal 
Ha intensity in a treatment that we call individual shell segmentation. The respective line ratios 
for the individual shells segmentation are given in tables |5] and [7] for IC 132 and in tables l9l and ITTI 
for BCLMP 93. The 12 + log(0/H) values, also for shell segmentation, are in figures l26landl34lfor 
IC 132 and BCLMP 93 respectively. 



4.2.6 About the abundance estimations 

The spatially resolved data on the M33 Hn regions allows the comparison of the different estima- 
tors on the assumption that the metal content of the ionized gas is well mixed giving a uniform 
total abundance across the Hn region. For IC 132 R 23 gives an average 12 + log(0/H) = 8.32 with 
standard deviation (s.d.)=0.27, the P method gives 12 + log(0/H) = 8.1 with s.d =0.30, the N2 
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Figure 15. 12 + log(CVH) abundance from the S23 indicator for IC 132 (left) and for the centre of M33 (right). 



method gives an average 12 + log(0/H) = 8.17 with s.d =0.27 and 03N2 gives 12 + log(0/H) = 
8.08 with s.d.=0.09. 

Hence, considering the mean of the abundance distribution we may say that the different indi- 
cators are consistent with each other in the sense that they give the same mean abundance value. 
Caution should arise for applying the P method to IC 132, given that the mean of the abundance 
distribution falls between the two branches and the turnaround region is not formally defined by 
the method. 

The results for the central Hn regions are equally in agreement for the mean abundance, but 
with a larger scatter in the three estimators. R 2 3 gives an average 12 + log(0/H) = 8.69 with 
s.d.=0.42, P gives 12 + log(0/H) = 8.43 with s.d.=0.38. The N2 method gives an average 12 + 
log(0/H) = 8.73 with s.d.=0.15 and 03N2 gives 12 + log(0/H) = 8.66 with s.d.=0.17. 

When the whole of the distribution is considered then the agreement between indicators seems 
to hold only for the brightest parts of the region and break badly for the spaxels at the borders, 
especially for the R23 and P methods. The uniform abundance assumption is strongly supported 
by 03N2 and N2, making these methods the preferred ones for abundance estimation for our data. 

© 2012 RAS, MNRAS 000,|2^?? 



32 Jesus Lopez-Hernandez et al. 

The P method for IC 132 is actually producing two separate distributions because the low and 
high abundance branches are selected for the same region. This is also the case for R23. If R23 is 
restricted to only one branch the abundance would be more uniform. This brings up the question 
about the validity of applying a branch selection for each spaxel, maybe in this case the selection 
by the [Nn]/[0n] criterion only applies to the integrated region. Another point to consider about 
empirical indicators is their sensitivity to the ionization parameter. From the maps it seems that 
03N2 gives low dispersion because it can effectively cancel out this effect. On the other hand R23 
shows high values especially in the lower left corner of IC 132 where we suspect that the diffuse 
ionized gas is dominant. The validity of applying R23 in such conditions is questioned. 

Similar results are found for the sulphur abundance. Using the S23 we find for IC 132 a sulphur 
abundance of 12 + log(S/H) = 6.79 with s.d.=0.28 and for the central Hn regions 12 + log(S/H) = 
7.02 with s.d.=0.34 

Regarding oxygen abundance, S23 also shows a relatively small scatter albeit higher than that 
of the oxygen abundance estimates using N2 and 03N2. 

When comparing direct abundance estimated with strong line methods, a nd using global mea- 



surements i.e. including the whole of the Hn regions. iKennicutt et al.1 (|2003r ) found that (0/H)^ 23 
is always higher than (0/H) re while the (O/H)^ method gives results in agreement with (0/H) re 
except for objects w ith high abundan ce values (12+log(0/H) T ' e > 8.5) where (O/H)^ is also 
larger than (0/H)r e - IStasihskal (I2002|) noted this systematic overestimates by strong line meth- 
ods, pointing out as possible causes the few number of objects with high quality T c determination 
and offsets between the temperatures assumed in the photoionization models and the temperatures 
observed. 

For N2, comparing the 



Denicolo et al. 



(2002) calibration with direct (O/H) determinations 



Perez-Montero & Diazl (|2005r) found that (O/H) n 2 overestimates abundances for low metallicity 
objects and underestimates them at high metallicity with a turnover at around 12 + log(0/H) Te ~ 
8.0. 

The selection of the empirical abundance indicator to use is primarily dictated by the available 
lines and their S/N, however in our case various indicators are available that produce considerable 
differences between the estimated abundance. The question of what empirical indicator can be 
trusted is not a simple one to answer. As the (O/H) calibrations are based on a set of objects 
with oxygen abundance derived with the T c method, the reliability of the relation depends on the 
quality, homogeneity and distribution of the reference objects. Ideally to construct an empirical 
estimator the reference sample should cover all the abundance range, in practice the reference is 
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generally weighted towards low or high abundance objects and with gaps at certain abundances. 
This situation is partially alleviated by the use of photoionization models, nevertheless its validity 
at different metallicities depends on our knowledge of transitional probabilities and collisional 
strengths for the different species. 

A warning may be raised about the use of empirical estimators for individual spaxels given that 
the methodology was developed for global observations. The comparison of the empirical abun- 
dance maps with the integrated zones shows in general good statistical agreement, and this may 
indicate that the estimators can be used for spatially resolved observations, some with more confi- 
dence (03N2,N2) than others (R23, P). Nevertheless determining the range of physical conditions 
where each empirical method is valid requires further investigation. 



5 DIAGNOSTIC DIAGRAMS IN 2D 

When interpreting emission-line spectra, it is important to be able to distinguish emission pro- 
duced by star-forming regions from other sources such as planetary nebulae, supernova remnants 
or even an active galactic nucleus. The conventional means for quantitatively classifying emission 
line objects and distingui shing between gas ionized by stars or by non-therm al processes are di- 



agnostic diagrams (BPT, iBaldwin et al.lll98ll : IVeilleux & Osterbrock Ill987|) involving the ratios 



of the strongest emission lines in the optical spectra. This traditional approach involves the com- 
parison of the integrated emission line ratios with global photoionization models leaving open the 
question of the validity of using these diagnostic diagrams with spatially resolved spectroscopic 
data. 

To search for systematic differences between the inner and outer Hll regions as well as detect- 
ing regions with line ratios indicative of non stellar ionization we computed the [Om]A 5007/H/?, 
[Nil] A 6584/Ha, [Sn]A 6717+6731/Ha, [On] A 3727/[0m]A 5007 and [On] A 3727/H/3 line ratios 
both spaxel by spaxel and the integrated values, to produce the distribution maps and investigate 
the emission line diagnostic diagrams and the excitation structure within the regions. 

Given the large ionization potential of [Olll] (54.9 eV) the [Om] A 5007/H/3 ratio traces highly 
ionized gas, and grows directly proportional to the degree of ionization. On the other hand because 
of the smaller ionization potential of both [Nil] and [Sn] (29.6 and 23.4 eV respectively) their line 
intensities trace the low ionization zones. These line ratio maps are shown in figure [[6]for IC 132 
and in figures [FT] and [T8] for the central regions. For IC 132 the "wall" in the lower left part is 
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present in most of the plots and is the strongest zone in the low ionization indicators maps ([Nil] 
A 6584/Ha, [Sll]A 67 17+673 1/Ha), pointing to the existence of diffuse ionized gas in this area. 

The Ha emission, being a tracer of star formation, can be expected to be spatially related to the 
position of the ionizing stars (this is true also for the highly ionized gas). Although the gas density 
and kinematics may play a role, in a first order approximation the maxima of Ha emission should 
follow the high ionization zones in the ionization ratio maps. From the [Olll]/H/3 ratio maps it is 
possible to see that this is true for the brightest part of IC 132 but not for the Hn regions in the 
central regions of M33. In particular the behaviour of the region BCLMP 93 is just the opposite, 
the brightest Ha peak coincides with a low ionization region and is partially surrounded by higher 
ionization gas. 
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Figure 18. Ionization sensitive ratios [OlI]A3727/[OlII]A5007 and [OlI]A3727/H/3 for the central region. 



Diagnostic BPT diagrams were created from the ionization and excitation sensitive line ratio 
maps shown before. In Fig.[19]the [Nil] A6584/H« vs [Olll] A5007/H/3 BPT diagrams are shown. 
The individual spaxels are plotted in the diagrams, coloured according to the Ha intensity, the 



insert shows t he map for reference. Bou ndary starburst model line s are shown as fo l 



dashed line is from 



Stasihska et al. 



line is that bv iKewlev & Do pita (2002), the dotted line is from Kauffmann et alJ (|2003h and the 



ows: the solid 



(2006). The location of the integrated spectra is indicated in 
black and white concentric circles. 

The main and obvious result from an inspection of Fig. [T9]is that there is an almost complete 
dichotomy in the distribution of points between the central and outer Hn regions. While most of 
the line ratios of the central region are vertically distributed on the right of the diagram, all of 
the points from IC 132 are located horizontally in the upper left part of the diagram. The regions 
occupied by the distributions have hardly any superposition, in fact, they are perpendicular to each 
other, but in both regions there are spaxels that clearly move away from the Hn region zone of the 
diagnostic diagram intruding in the AGN zone, mostly the Seyfert zone in the case of IC 132 and 
the LINER and transition zone for the M33 central region. Curiously, no anomalous line ratios are 
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detected for the position of individual spaxels near the M33-X8 X-ray source in the central region. 
The position of the emission line zone nearest to the X-ray source is indicated with a triangle 
symbol in the left panel of figure [191 
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-2.0 -1.5 -1.0 -0.5 -2.0 -1.5 -1.0 -0.5 

log([Nll]A6584/Ha) log([Nll]A6584/Ha) 

Figur e 19. BPT diagnostic diagrams for each spaxel in the exte rnal zone IC 132 (left) and the centra l zone s of M33 (right). Different boundaries for photoionization are represented: continuous line from 
iKewlev & Dopital <2002h . dotted line from lKauffmann etalj 120031) and dashed line from IStasihska et"aD feOOfj) . The colours of individual squares indicate the spatial position according to the inserted map of 
the regions. Black and white concentric circles are the position of the integrated field for IC 132 and BCLMP 93. Symbols connected by coloured lines are the results from photoionization models computed at 
increasing distance from the ionization source, where the symbols closest to the boundary lines are the final result of the computation. Three temperatures are given for a combination of log U and abundance. For IC 
132 the best fit was found with models running horizontally in the diagram (green lines), while the central zone was fitted not only by horizontal but specially by vertical (orange lines) models. The triangle symbol 
in the diagram for the centre indicates the emission line zone closer to the strong X-ray source (see text). 
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5.1 Matching the diagnostic diagrams to pseudo-3D photoionization models 

Bearing in mind that we are dealing with spatially resolved data and in an effort to understand the 
dichotomy in the distribution of individual spaxel line ratios in the BPT diagrams shown in Fig.[T9l 
we have computed what can be called pseudo-3D photoionization models. Our grid of models 
have spherical closed geometry, constant density with a black body as the ionizing source and a 
range of ionization parameters (log U) from -2.0 to 0.0, metallicities (Z) from 0.1 to 2.OZ and 
stellar temperature (T*) from 30,000 to 50,000 K. As a first approach to full 3D models we have 



used all intermediate shell results from each Cloudy (|Ferland et al] |1998) photoionization model 
to represent the radial behaviour of the relevant gas parameters inside the nebula.The models are 
discussed in detail in a forthcoming paper (Lopez et al. in preparation). 

Fig. [H] shows in green-black solid lines the results from the photoionization models. Each line 
corresponds to a model with a given combination of log U, Z and T*, following the computed 
results at increasing distance from the inner face of the cloud. The end or outer edge of the com- 
putation is marked with the corresponding symbol in bold and is generally the one closest to the 
boundary lines. 

It is possible to see in Fig. [19] that the models are capable of reproducing the position and 
spread of values of both regions even though the spaxels in IC 132 show a horizontal dispersion 
while those of the central zones have a vertical dispersion. In particular, the models that best fit the 
data for IC 132 have log(U)=-l, T> 35,000 and 40,000°K and Z=O.1Z . T> 35,000 fits the high 
brightness spaxels and 1^=40,000 lies above these spaxels but fits better the low brightness ones. 
T*=30,000 is below most of the spaxels but fits a small spur of objects around [Nll]/Ha=-1.05. 

For the central zone, models with the same temperatures (3.0, 3.5 and 4.0 kK), log(U)=-l and 
Z=2.0Z Q span the location of most of the individual spaxels. Another set of models with the same 
T*, but log(U)=-2 and Z=0.1Z o also fit the observations but with vertical spread; these are the 
orange broken lines in the right panel of figure [19] These vertical models only fit the rightmost 
distribution of central spaxels, corresponding to the faint ones. The bright spaxels fall on the lower 
left tip of the distribution cloud and could not be fitted by our set of models. However the integrated 
position for BCLMP 93 (concentric circles in the figure) is well fitted by both sets of models. 

For the same temperatures, the position of the horizontal models is controlled by the metallic- 
ity, the higher it is the lower the [Om]5007/H/5 ratio. The orientation of the models is controlled 
by the ionization parameter, being vertical for low values of U. 

For the central zone we favour the vertical models as indicators of the physical conditions in 
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the nebula, mainly because they are consistent with the high metallicity expected in this region, 
and the horizontal models instead would represent regions of low metallicity. 

As mentioned above while the spaxels with high Ha surface brightness are inside the boundary 
starburst line of the BPT diagram, many of the low Ha surface brightness spaxels that are further 
away from the ionizing source appear outside the boundary; this obviously doesn't mean that an 
AGN exists. Possible explanations are that those pixels have low S/N or perhaps diffuse radiation 
has to be considered. In any case, this raises a question about the wisdom of using photoionization 
models computed for integrated systems to define a reference boundary for photoionization and 
apply it to diagnostic diagrams for resolved systems. 



6 WOLF-RAYET STARS 



Wolf-Rayet (WR) stars are hot and lumin ous, with broad emissi on lines in the optical range due to 
a considerable mass loss through winds (|Abbott & Contilll987|) . They represent an evolved stage 
of O stars when the radiation pressure can not be counterbalanced by gravity and the outer layers of 
the stars are blown away. WR more common subtypes are the WN that show a prominent nitrogen 
emission in the 4600-4720 A range, known as the blue bump (BB), and the WC, with carbon 
emission in the 5750-5870 A range, known as the red bump (RB). The [Hell]A4686A emission 
line is present in both types. 

If Wolf-Rayet features are detected then stars m ore massive than M wr ~ 25M© must be 



present in the cluster if solar metallicity is presumed (IPindao et al. 



2002|) . Assuming an instanta- 



neous burst, the WR features are visible f rom roughly 3 Myr to a lmost 7 Mys of age, their lifetimes 



getting shorter with lower metallicity (jSchaerer & Vacca 



1998). The WR feature can serve as a 



co nstraint to estimate the age of th e burst if its metallicity is known. 



D'Odorico & Benvenuti 



(119831) studied the integrated UV (IUE) and optical (ESO 3.6m) spec- 
trum of IC 132. They found the presence of the blue bump (A 4686), indicating the existence of 
WR stars, which they classified as WN4. They estimated the ionizing source as a single star with 
T,=40,000±5,000°K and mass larger than 100M© or a system of 3 components (04V+09.5I+WN4). 
They noted the need for spatial resolution to determine the extent and nature of the WR stars (few 
over luminous objects vs. many normal WR). 

We used the IFS observations to map the spatial distribution of the WR features by measuring 
the detected BB in each individual spectrum and reconstructing the spatial map in the same fash- 
ion as the emission line maps were constructed. We detected two peaks in the BB distribution in 
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Table 3. Blue bump luminosity, EW and BB/H/3 for the two regions identified in IC 132. 



region 


log(L) 


EW(BB) 


BB/H/3 




erg s — 1 


A 




A 


37.84 


15 


0.38 


B 


37.89 


13 


0.41 



IC 132. Fig. [20] shows the map of the blue feature flux measured by fitting a gaussian in a 70A 
interval centred at 4686A, with the Ha contour overlaid. Also shown is a zoom of the BB emission 
in the region where we detect two peaks and the integrated spectra of both regions. The BB lu- 
minosity, equivalent width and fl ux normalized to H /3 are shown in table |3] Assuming an average 



luminosity for the BB of 10 36 5 (|Pindao et al 



2002L for a single WN7 star) the estimated number 



of WN stars in region A is 21 and in region B is 24. Taking the total number of Lyman photons 
estimated bv lD'Odorico & Benvenutil (11983) N /-,,, = 8.2 x 10 49 s _1 and assuming that a single WR 



produces 1.7 x 10 49 s 1 (|Schaerer & Vacca 



1998 



WN), the derived number of stars is 5 for the 



whole region. It is interesting to note the discrepancy found in the number of WR computed by 
different methods. This discrepancy may be linked to aperture effects in the observations. 
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Given that the models of WR formation predict that as a result of the increased mass loss rates 



the n umber of WR will increase with increasing metallicity (|Schaerer & Vacca 



1998; M evnet & Maeder 



20051), it is surprising that no WR features are detected in the higher metallicity central region of 



M 33. One has to consider that the spectrum of the central region BCLMP 93 is of lower S/N 
than that of IC 132 which makes more difficult the detection and measurement of the WR bumps. 
Nevertheless, a limit to the number of WR stars in BCLMP 93 that may be present but can not 
be confidently detected can be obtained by estimating the blue bump with different continuum 
placements. In this way we find that BCLMP 93 could contain at most 4 WN7 stars. 



7 SEGMENTATION ANALYSIS 

Although the intensity maps provide a view of the spatial distribution of the emission lines flux, 
temperatures, abundances and other estimated values, an easier quantitative representation of vari- 
ations with position as well as a study of the effect of aperture size is provided by analyzing the 
Hll region in concentric layers. To this end we have adopted an approach that defines segments 
or shells of the Hll region according to the Ha surface brightness level. The segmentation of the 
observed FOV allows the study of the faint areas where the spaxel by spaxel analysis can not be 
performed given the low S/N of individual elements. However when the flux is integrated in a 
single spectrum the S/N improves and conditions can be estimated, albeit with the loss of spatial 
azimuthal resolution. This kind of analysis is a clear advantage of IFS over long slit observations 
for the study of faint areas. 

We present the results in two forms: 

I - The individual shell values, i.e. considering each section as an individual region, and 

II - The cumulative values, i.e. adding the sections in a cumulative way from the centre to the 
outer boundary, going from AB=A+B to AJ=A+B+ . . . +I+J. 

Therefore in this notation A is the spectrum of the brightest central spaxels, while AJ repre- 
sents the integrated spectrum of the whole region (see figure [2T| for the segmentation of IC 132). 
The same notation is used for BCLMP 93 but with lowercase letters going from a to h (see figure 

7.1 Integrated properties of IC 132 

To analyze the radial properties of IC 132 the field was segmented in 10 sections of equal Ha 
surface brightness covering from the brightest to the faintest spaxels. The shells were labeled from 
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A RA [pix] 

20 40 60 




20 40 60 




Figure 21. IC 132 Ha map with concentric apertures. The single letter (A through J) designates the shells area. 

A to J as shown in Fig. [2U The segment sizes were defined to have constant integrated flux in Ha 
after applying the extinction correction. 

As can be seen from tables 4 and 8 the integrated flux of H/? is roughly the same for all sections 
in both regions, i.e. about 80x 10~ 13 erg s _1 cm~ 2 for IC 132 and about 20x 10~ 13 erg s _1 cm~ 2 
for BCLMP 93; small variations are present because of the discrete jumps in flux between spaxels 
and no interpolation was used to define exactly the isocontours. 

The emission line measurements, their errors, and the resulting line ratios, temperatures and 
abundances, are given in tables 4 to 7 for the individual and cumulative spectra respectively. 
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Table 4. Reddening corrected integrated fluxes forIC 132 individual shells. Fluxes are normalized to H/3 = 100 



4^ 



Individual 



Ion 



3727 


[Oil] 


150.7±8.0 


164.7±8.9 


192.4±11.7 


210.5±15.7 


204.5±16.5 


188.9±17.6 


221.6±21.0 


241.3±27.7 


267.2±44.8 


634.6±282.6 


4010 


H<5 


33.5±2.4 


31.5±2.6 


31.9±3.1 


32.0±3.9 


31.5±3.9 


31.0±4.0 


30.5±4.8 


29.5±7.0 


28.4±10.8 


24.6±54.3 


4340 


H7 


JU.U±2.8 


48.5±2.6 


4y.z±j.o 


49.3±5.3 


48.2±5.6 


47.7±6.7 


47. 1 ±7.8 


4j./±lU.U 


43. 0± 15. 1 


31.4±52.5 


4363 


[O III] 1 


6.7±1.2 


6.4±1.3 


7.7±1.7 


8.4±2.3 


10.7±2.7 


11.2±3.3 


14.0±4.3 


18.9±6.3 


29.0±11.6 


127.0±67.6 


4861 


tm 


100.0±3.1 


100.0±4.1 


100.0±5.7 


100.0±7.9 


100.0±7.9 


100.0±8.6 


100.0±9.6 


100.0±12.1 


100.0±18.5 


100.0±48.2 


4959 


[O III] 


175.1±6.4 


171.4±7.3 


164.7±8.0 


159.6±11.8 


157.8±11.6 


158.0±12.3 


154.7±13.4 


155.4±17.4 


160.9±27.5 


208.3±85.4 


5007 


[O III] 


500.1±12.2 


485.6±15.4 


463.5±20.3 


449.1 ±27.8 


443.5±28.1 


445.1 ±30.8 


435.0±33.3 


433.3±42.4 


444.6±66.5 


576.1 ±209.9 


6312 


[S III] 


2.7±1.3 


3.4±2.1 


3.9±2.7 


4.6±3.4 


4.9±4.0 


5.5±4.7 


6.1±5.8 


7.4±7.9 


12.0±13.8 


51.9±60.3 


6548 


[Nil] 


2.9±1.1 


2.8±1.2 


3.2±1.3 


3.4±1.8 


3.5±2.2 


3.5±2.1 


4.1±2.7 


4.7±3.4 


6.7±4.4 


16.4±14.8 


6564 


Ha 


284.5±7.1 


284.1±8.8 


284.3±12.2 


284.2±17.8 


284.8±18.8 


284.5±19.5 


284.5±22.5 


284.8±28.0 


285.8±40.0 


317.8±112.4 


6584 


[Nil] 


6.7±0.8 


7.7±0.9 


9.1±1.2 


9.8±2.2 


9.6±2.7 


9.1±2.4 


10.1±3.3 


11.3±4.2 


12.7±5.0 


20.0±15.7 


6717 


[S II] 


8.9±0.8 


9.9±0.9 


11.7±1.2 


12.7±3.2 


12.3±4.2 


11.5±3.0 


12.9±4.6 


14.0±5.4 


16.1±4.9 


26.0±16.2 


6731 


rsm 


5.7±0.4 


6.6±0.7 


7.6±0.9 


8.4±1.5 


8.2±1.7 


7.3±2.0 


8.6±2.5 


9.2±3.4 


9.5± 5.5 


11.6± 18.8 


7136 


[Ar III] 


8.3±1.0 


6.6±1.2 


7.7±1.4 


8.3±2.9 


8.3±1.8 


8.4±2.2 


8.2±3.0 


7.3±4.8 


5.6±4.3 


2.9±12.6 


9068 


[Sill] 


27.5±6.3 


24.0±5.4 


28.4±6.5 


30.8±8.2 


30.6±7.6 


28.1 ±7.5 


28.2±9.3 


23.4±9.9 


18.9±13.2 


13.1±33.5 


I(H/3)[erg s^cm" 2 10 13 J 


79.29±0.18 


82.86±0.24 


79.24±0.32 


80.21 ±0.45 


80.99±0.45 


80.17±0.49 


80.52±0.55 


80.50±0.69 


80.21±1.05 


71.20±2.43 






0.50±0.06 


0.42±0.04 


0.40±0.04 


0.40±0.05 


0.37±0.05 


0.36±0.05 


0.36±0.06 


0.39±0.06 


0.49±0.09 


1.30±0.13 


EW(H/3)[A] 


130.34± 11.42 


139.98±11.83 


194.10±13.93 


292.21±17.09 


284.06±16.85 


266.77±16.33 


218.40±14.78 


183.85±13.56 


132.50± 11.51 


23.46±4.84 



PS 

I 

TO 
TO 

3 
s 

TO 

<N 

TO 
«-+ 

& 



1 Contaminated by mercury street lamps. 




* -0 

* -0 



Table 5. Integrated temperatures, abundances and line ratios for IC 132 individual shells. 



Individual 





A 


B 


C 


D 


E 


F 


G 


H 


T e ([0lll]) a 


12800±900 


12700±1000 


13900±1300 


14600±1800 


16500±2100 


16800±2600 


19200±3500 


23400±5700 


Te([Oll])pM03 


11300±400 


12100±300 


12800±400 


13200±500 


14200±500 


14400±600 


15600±700 


17300±900 


T e ([Oll]) G92 


12000±600 


11900±700 


12700±900 


13200±1300 


14500±1500 


14800±1800 


16500±2400 


19400±4000 


T e ([SlIl]) 63 12 


10700±100 


12800±200 


12500±300 


13100±400 


13500±500 


15100±700 


16500±1000 


22400±2800 


Te([SlIl]) S o 23 


9300±1900 


9600±1900 


9200±1700 


9000±1500 


900±1500 


9100±1600 


9100±1600 


9600±1700 


re([Sni]H08 


12000±1000 


11900±1200 


13300±1600 


14000±2200 


16400±2500 


16800±3100 


19800±4100 


24700±6800 


12+log(0+/H+)<> 


7.46±0.19 


7.50±0.22 


7.46±0.26 


7.44±0.33 


7.29±0.32 


7.23±0.39 


7.16±0.42 


7.01±0.52 


12+log(0++/H+) b 


7.94±0.19 


7.93±0.23 


7.80±0.26 


7.73±0.32 


7.60±0.30 


7.58±0.36 


7.44±0.38 


7.27±0.45 


12+log(0/H) b 


8.06±0.27 


8.07±0.32 


7.97±0.37 


7.91±0.46 


7.77±0.44 


7.74±0.53 


7.62±0.57 


7.46±0.68 


12+log(S+/H+) 


5.49±0.32 


5.35±0.36 


5.46±0.44 


5.46±0.56 


5.43±0.66 


5.29±0.79 


5.29±1.07 


5.14±1.73 


12+log(S++/H+) 


6.69±0.18 


6.49±0.17 


6.58±0.18 


6.57±0.22 


6.55±0.19 


6.45±0.19 


6.39±0.26 


6.12±0.35 


ICF(S + + S^ + ) 


1 \ f\ AO 

l.3l±0.08 


1.28±0.10 


l.22±0.ll 


1.19±0.13 


1.20±0.13 


l.22±0.l5 


1.19±0.16 


I.l8±0.l8 


12+log(S/H) 


6.83±0.37 


6.63±0.39 


6.70±0.47 


6.68±0.60 


6.66±0.68 


6.56±0.8l 


6.49±1.10 


6.24±l.77 




-l 34±0 02 


-1 55±0 02 


-l 35±0 02 


-1 31+0 03 


-1 19±0 03 


-l 27±0 03 


-1 20±0 03 


-l 29±0 04 


12+log(N+/H+) 


5 92±0 01 

.. • • y — 1 v / ■ v./ 1 


5 98±0 01 


5 99±0 01 


5 98±0 02 

~ ' • / V l 1 v./ . \ / — 


5 88±0 03 


5 85±0 03 


5 80±0 03 


5 73±0 04 




-1 54±0 06 


-1 52±0 06 

J- ■ -y 1 V / . V / \ I 


-1 47±0 06 

-L • T / 1 \J ■ \J\J 


-1.46±0.07 


-1.41±0.07 


-l 38±0 08 


-1 36±0 08 


-1 28±0 09 


12+log(0/H) (N2) 


8.02±0.06 


8.05±0.09 


8.08±0.20 


8.09±0.43 


8.09±0.48 


8.08±0.50 


8.10±0.71 


8.12±1.19 


12+l0g(0/H) Slower) 


8.06±0.08 


8.07±0.10 


8.10±0.13 


8.11±0.18 


8.10±0.18 


8.07±0.20 


8.12±0.22 


8.15±0.27 


12+log(0/H) (Ar 3 3 ) 


8.03±0.06 


7.91±0.09 


8.03±0.08 


8.09±0.14 


8.10±0.09 


8.10±0.11 


8.10±0.15 


8.04±0.29 


12+log(0/H) (S3O3) 


8.27±0.05 


8.24±0.05 


8.29±0.04 


8.32±0.05 


8.32±0.04 


8.31 ±0.05 


8.31 ±0.06 


8.27±0.09 


RS2 


1.21±0.03 


1.45±0.03 


1.49±0.03 


1.51±0.05 


1.71±0.06 


1.62±0.05 


1.86±0.06 


2.21±0.07 


log(S0 23 ) 


-0.77±0.26 


-0.82±0.26 


-0.74±0.25 


-0.71 ±0.24 


-0.71 ±0.24 


-0.73±0.24 


-0.74±0.23 


-0.81±0.23 


N2 


-1.63±0.01 


-1.57±0.02 


-1.49±0.05 


-1.46±0.11 


-1.47±0.12 


-1.50±0.12 


-1.45±0.18 


-1.40±0.31 


log(R23) 


0.92±0.04 


0.91±0.05 


0.91±0.06 


0.91±0.09 


0.91±0.09 


0.90±0.10 


0.91±0.11 


0.92±0.14 




0.65±0.06 


0.60±0.06 


0.51±0.07 


0.46±0.09 


0.47±0.10 


0.50±0.11 


0.43±0.11 


0.39±0.14 


log(Ar 3 3 ) 


-1.78±0.04 


-1.87±0.06 


-1.78±0.06 


-1.73±0.12 


-1.73±0.08 


-1.73±0.10 


-1.72±0.13 


-1.77±0.23 


log(S 3 3 ) 


-1.15±0.08 


-1.19±0.07 


-1.10±0.08 


-1.06±0.09 


-1.05±0.08 


-1.08±0.08 


-1.07±0.11 


-1.15±0.15 


log(u) ([SiiMSm]) 


-1.46±0.16 


-1.58±0.15 


-1.61±0.16 


-1.62±0.21 


-1.59±0.21 


-1.55±0.21 


-1.62±0.26 


-1.81±0.31 


log(u) ([OiiMOm]) 


-2.50±0.02 


-2.54±0.02 


-2.61±0.02 


-2.65±0.03 


-2.65±0.03 


-2.62±0.04 


-2.68±0.04 


-2.71±0.05 


log(r?') 


0.26±0.07 


0.24±0.06 


0.31±0.07 


0.35±0.08 


0.36±0.09 


0.35±0.10 


0.39±0.10 


0.32±0.12 


log([Oin]A5007/H/3) 


0.70±0.01 


0.69±0.01 


0.67±0.02 


0.65±0.03 


0.65±0.03 


0.65±0.03 


0.64±0.03 


0.64±0.04 


log([Nn]A6584/Ha) 


-1.63±0.05 


-1.57±0.05 


-1.49±0.06 


-1.46±0.09 


-1.47±0.12 


-1.50±0.11 


-1.45±0.14 


-1.40±0.16 


log([Sn]AA6717,31/Ha) 


-1.27±0.06 


-1.25±0.06 


-1.17±0.06 


-1.13±0.08 


-1.15±0.10 


-1.20±0.09 


-1.16±0.10 


-1.13±0.10 



s 
s 



4^ 



a Upper limit. Lower limit. 



Table 6. Reddening corrected fluxes for IC132 accumulated shells. Fluxes are normalized to H/3 = 100 



Cumulative 



A 


Ion 


A A 

A A 


A Li 


A I 


A TO 
AD 


AE 


Ar 


ALr 


A I r 
All 


Al 


A.J 


3727 


[0 II] 


1 jU. /±o.U 


I J / .o±o.U 


ioy.z±o. / 


i /y.j±y.j 


1 o4.J± 111. 1 


loj.Z±lU.O 


1 on A J— 1 1 A 

iyu.4± 1 1 .4 


lyo. /± i j.u 


7H/1 1 < < 
ZU4. J± 1 J.J 


7/17 Q^77 Q 
Z4Z.O±JJ.y 


4010 


H() 


33.5±2.4 


32.5±2.4 


3z.j±z.o 


77 1_L1 

3z.z±Z.o 


77 i j_7 n 


31.y±3.0 


7 1 7X1 7 


31.4±3.7 


31.1i4.5 


30.o±y.4 


4340 


H 7 


50.0±2.8 


49.2±2.2 


49.2±2.5 


49.2±2.9 


49.0±3.3 


48.8±3.8 


48.5±4.3 


48.2±5.0 


47.6±6.2 


45.7±11.4 


4363 


[0 III] 1 


6.7±1.2 


6.6±1.0 


6.9±1.1 


7.3±1.3 


8.0±1.5 


8.5±1.8 


9.2±2.1 


10.4±2.5 


12.4±3.3 


22.6±7.8 


4861 


[H/3] 


100.0±3.1 


100.0±3.5 


100.0±4.0 


100.0±4.7 


100.0±5.2 


100.0±5.6 


100.0±6.1 


100.0±6.8 


100.0±8.2 


100.0±12.8 


4959 


[O III] 


175.1±6.4 


173.2±6.7 


170.4±6.8 


167.7±7.6 


165.7±8.2 


164.4±8.7 


163.0±9.2 


162.0±10.1 


161.9±12.0 


166.1±19.3 


5007 


[O III] 


500.1±12.2 


492.7±13.4 


483.1±15.1 


474.6±17.3 


468.3±19.0 


464.5±20.6 


460.2±22.2 


456.8±24.7 


455.4±29.8 


466.4±47.9 


6312 


[S III] 


2.7±1.3 


3.1±1.7 


3.4±2.0 


3.7±2.4 


3.9±2.7 


4.2±3.0 


4.4±3.4 


4.8±4.0 


5.6±5.1 


9.7±9.8 


6548 


[N II] 


2.9±1.1 


2.8±1.1 


3.0±1.2 


3.1±1.3 


3.2±1.4 


3.2±1.5 


3.3±1.7 


3.5±1.9 


3.8±2.1 


5.0±3.1 


6564 


Ha 


284.5±7.1 


284.3±7.7 


284.3±8.8 


284.3±10.4 


284.4±11.7 


284.4±12.7 


284.4±14.0 


284.4±15.7 


284.6±18.4 


287.8±28.2 


6584 


[N II] 


6.7±0.8 


7.2±0.8 


7.8±0.9 


8.3±1.1 


8.6±1.4 


8.7±1.6 


8.9±1.8 


9.2±2.1 


9.6±2.4 


10.5±3.4 


6717 


[S II] 


8.6±1.9 


9.2±1.4 


10.0±0.9 


10.7±1.1 


ll.Oil.O 


11.3±2.2 


11.6±1.1 


11.7±1.0 


12.2±2.8 


12.8±1.5 


6731 


[S II] 


5.8±0.3 


6.3±0.3 


6.2±0.5 


6.3±0.6 


6.6±0.8 


6.6±0.9 


6.9±1.0 


6.9±1.2 


7.2±1.4 


7.0±2.0 


7136 


[Ar III] 


8.3±1.0 


7.4±0.9 


7.5±1.0 


7.7±1.3 


7.8±1.3 


7.9±1.3 


7.9±1.3 


7.9±1.5 


7.6±1.7 


7.2±2.1 


9068 


[Shi] 


27.5±6.3 


25.7±5.0 


26.6±5.3 


27.6±5.8 


28.2±5.9 


28.2±5.9 


28.2±6.3 


27.6±6.5 


26.6±6.9 


25.1±8.9 


UHPftergs- 1 ™- 2 10~ 13 ] 


79.29±0.18 


162.15±0.40 


241.39±0.69 


321.61±1.07 


402.61±1.47 


482.78±1.91 


563.35±2.43 


643.90±3.10 


724.18±4.18 


795.15±7.19 


C(H/3) 




0.50±0.06 


0.45±0.03 


0.43±0.03 


0.42±0.02 


0.40±0.02 


0.40±0.02 


0.39±0.02 


0.39±0.02 


0.40±0.02 


0.49±0.02 


EW(H/3)[A] 


130.34±11.42 


135.09±11.62 


150.07±12.25 


170.79±13.07 


185.69±13.63 


195.56±13.98 


198.54±14.09 


196.59±14.02 


186.62± 13.66 


114.99± 10.72 



Contaminated by mercury street lamps. 



Table 7. Integrated temperatures, abundances and line ratios for IC 132 accumulated shells. 



Cumulative 


































AA 


AB 


AC 


AD 


AE 


AF 


AG 


AH 


AI 


AJ 


Te([Om]) 


12700±900 


12700±700 


13100±900 


13400±1000 


14000±1200 


14500±1300 


15000±1600 


16000±1900 


17600±2600 


25300±6960 


Te([Oll])pM03 


12100±300 


12100±200 


12300±200 


12600±300 


12900±300 


13100±300 


13500±400 


14000±400 


14800±500 


18000±1000 


r e ([On]) G92 


11900±600 


11900±500 


12200±600 


12400±700 


12800±800 


13100±900 


13600±1100 


14200±1300 


153O0±1800 


20700±4900 


r e ([SlIl]) 63 l2 


10600±100 


11600±200 


12000±200 


12200±200 


12500±200 


13O00±30O 


13400±400 


14300±500 


16200±800 


27900±4500 


r e ([Sm]) sc , 23 


9400±1900 


9500±19OO 


9400±1800 


9300±1800 


9200±1700 


9100±1700 


9200±1700 


9200±1700 


9300±1700 


9600±1700 


TpfrSlIllr/nK 


12000±1000 


12000±900 


124000±1000 


12800±1200 


13500±1400 


14000±1600 


14700±1900 


15900±2200 


17700±3000 


27000±8300 


12+log(0+/H+) 


7.45±0.19 


7.47±0.16 


7.47±0.18 


7.46±0.21 


7.43±0.22 


7.39±0.25 


7.35±0.27 


7.30±0.30 


7.22±0.36 


6.94±0.56 


12+log(0++/H+) 


7.94±0.19 


7.94±0.17 


7.89±0.19 


7.85±0.21 


7.80±0.22 


7.76±0.24 


7.70±0.26 


7.64±0.28 


7.54±0.33 


7.24±0.48 


12+log(0/H) 


8.06±0.27 


8.06±0.23 


8.03±0.26 


8.00±0.30 


7.95±0.31 


7.91±0.34 


7.87±0.37 


7.80±0.41 


7.71±0.48 


7.41±0.74 


12+log(S+/H+) 


5.49±0.32 


5.43±0.29 


5.44±0.32 


5.45±0.37 


5.44±0.41 


5.42±0.45 


5.40±0.53 


5.35±0.64 


5.27±0.88 


4.96±1.92 


12+log(S++/H+) 


6.69±0.18 


6.59±0.14 


6.58±0.15 


6.58±0.16 


6.57±0.15 


6.55±0.15 


6.52±0.16 


6.47±0.17 


6.37±0.20 


6.04±0.27 


ICF(S+ + S++) 


1.31 ±0.08 


1.29±0.07 


1.27±0.08 


1.25±0.09 


1.24±0.09 


1.23±0.10 


1.22±0.11 


1.22±0.12 


1.21±0.14 


1.20±0.20 


12+log(S/H) 


6.84±0.37 


6.73±0.32 


6.72±0.35 


6.71 ±0.40 


6.69±0.43 


6.67±0.48 


6.64±0.55 


6.59±0.67 


6.48±0.90 


6.15±1.93 


log(S/0) 


-1.34±0.02 


-1.45±0.02 


-1.42±0.02 


-1.39±0.02 


-1.35±0.02 


-1.33±0.02 


-1.31 ±0.02 


-1.30±0.03 


-1.31 ±0.03 


-1.34±0.04 


12+log(N+/H+) 


5.92±0.01 


5.95±0.01 


5.96±0.01 


5.97±0.01 


5.95±0.02 


5.93±0.02 


5.91±0.02 


5.89±0.02 


5.84±0.02 


5.65±0.03 


log(N/0) 


-1.54±0.06 


-1.53±0.05 


-1.51 ±0.05 


-1.49±0.05 


-1.47±0.05 


-1.46±0.06 


-1.44±0.06 


-1.42±0.07 


-1.38±0.07 


-1.29±0.10 


12+log(0/H) (N2) 


8.02±0.06 


8.04±0.07 


8.05±0.09 


8.07±0.13 


8.07±0.17 


8.07±0.21 


8.08±0.25 


8.08±0.32 


8.09±0.46 


8.11±1.14 


12+l0g(0/H) (R23Wer) 


8.06±0.08 


8.07 ±0.08 


8.08±0.10 


8.09±0.11 


8.09±0.12 


8.09±0.13 


8.09±0.14 


8.1O±0.16 


8.11±0.19 


8.19±0.30 


12+log(0/H) (Ar 3 3 ) 


8.03±0.06 


7.98±0.06 


8.00±0.06 


8.02±0.08 


8.03±0.08 


8.04±0.07 


8.05±0.08 


8.05±0.09 


8.04±0.10 


7.99±0.14 


12+loa(0/H) (S iOO 


8.27±0.05 


8.25±0.04 


8.27±0.04 


8.28±0.04 


8.29±0.04 


8.29±0.04 


8.30±0.04 


8.29±0.04 


8.28±0.05 


8.26±0.07 


RS2 


1.50±0.03 


1.47±0.02 


1.61±0.01 


1.69±0.02 


1.66±0.01 


1.73±0.03 


1.68±0.02 


1.70±0.02 


1.71 ±0.04 


1.83±0.02 


log(S0 23 ) 


-0.78±0.26 


-0.80±0.26 


-0.78±0.25 


-0.76±0.25 


-0.75±0.25 


-0.75±0.25 


-0.74±0.25 


-0.75±0.24 


-0.77±0.24 


-0.81±0.23 


N2 


-1.63±0.01 


-1.60±0.01 


-1.56±0.02 


-1.53±0.03 


1 CO 1 f\ f \ A 

-1.52±0.04 


-1.52±0.05 


-1.51 ±0.06 


-1.49±0.08 


-1.47±0.1 1 


i a a i a 

-1.44±0.29 


log(R 23 ) 


0.92±0.04 


0.92±0.04 


0.92±0.05 


0.91±0.05 


0.91 ±0.06 


0.91 ±0.06 


0.91 ±0.07 


0.91±0.08 


0.91 ±0.09 


0.94±0.15 


l°g(0 3 2) 


0.65±0.06 


0.63±0.06 


0.59±0.06 


0.55±0.06 


0.54±0.06 


0.53±0.07 


0.52±0.07 


0.50±0.08 


0.48±0.09 


0.42±0.16 


log(Ar 3 3 ) 


-1.78±0.04 


-1.82±0.04 


-1.81±0.05 


-1.79±0.06 


-1.78±0.06 


-1.77±0.06 


-1.76±0.06 


-1.76±0.07 


-1.78±0.08 


-1.81±0.11 


log(S 3 3 ) 


-1.15±0.08 


-1.17±0.06 


-1.15±0.06 


-1.12±0.07 


-1.11±0.07 


-1.10±0.07 


-1.10±0.07 


-1.11±0.08 


-1.12±0.09 


-1.15±0.12 


log(u) ([Sll]/[SIII]) 


-1.42±0.17 


-1.52±0.12 


-1.53±0.12 


-1.54±0.13 


-1.55±0.12 


-1.56±0.15 


-1.58±0.13 


-1.60±0.14 


-1.66±0.19 


-1.70±0.22 


log(u) ([Oll]/[Olll]) 


-2.50±0.02 


-2.52±0.02 


-2.55±0.02 


-2.58±0.02 


-2.59±0.02 


-2.60±0.02 


-2.61 ±0.02 


-2.62±0.03 


-2.64±0.03 


-2.69±0.06 


log(V) 


0.28±0.07 


0.25±0.06 


0.28±0.06 


0.31 ±0.06 


0.32±0.06 


0.32±0.06 


0.32±0.07 


0.33±0.07 


0.31±0.08 


0.35±0.14 


log([Om]A5007/H/3) 


0.70±0.01 


0.69±0.01 


0.68±0.01 


0.68±0.02 


0.67±0.02 


0.67±0.02 


0.66±0.02 


0.66±0.02 


0.66±0.03 


0.67±0.04 


log([Nll]A6584/Ha) 


-1.63±0.05 


-1.60±0.05 


-1.56±0.05 


-1.53±0.06 


-1.52±0.07 


-1.52±0.08 


-1.51±0.09 


-1.49±0.10 


-1.47±0.10 


-1.44±0.14 


log([SlI]AA6717,31/Ha) 


-1.30±0.06 


-1.27±0.04 


-1.24±0.03 


-1.22±0.03 


-1.21±0.03 


-1.20±0.06 


-1.19±0.03 


-1.18±0.03 


-1.17±0.07 


-1.16±0.05 



References: PM03 lPerez-Montero & Dfazl <2003h : G92 lGarnetJ jl992l) : H06 lHagele et alj J2006h 
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The resulting integrated properties for IC 132 are shown in plots [22] to [30] for electron temper- 
atures, ionic and total abundances, empirical abundances, line ratios, ionization sensitive ratios, 
extinction and EW(H/?). Each figure shows the values for individual shells on the left and for cu- 
mulative segmentation on the right panels. Individual shells I an J are not shown given the large 
errors obtained. Errors were calculated by formal propagation of each line flux associated error. 
For T e ([Oll])pAf03, also the error in n c determination at the low density limit was added in quadra- 
ture, to take into account the spread of the models in the original paper. 

Figure [22] shows the radial trend of the electron temperatures. Inside the errors there is a clear 
tendency for the temperature to increase with radius, except for T e ([Slll]) computed with the S0 2 3 
method, which shows no temperature gradient. The T c increasing tendency is reflected in the radial 
trend of abundance shown in Fig. [23] where it is possible to see that the estimated O/H goes 
to smaller abundances in the outer regions. However, given the large errors obtained, the radial 
behaviour is also compatible with a constant value for the whole region. In such case the total 
(O/H) may be set between 7.8 and 8.2. 

For figure [24] although a slight decrease with radius is noted, the ionic and total abundances 
are compatible with a constant value within errors. Total (S/H) is between 6.4 and 7.1. The mean 
value from the 2D map of 6.62 may as well be adopted as the representative value of the region. 

Interestingly, in figure [25] S/O has a rather flat trend while N/O displays an apparent increase 
with radial distance in the individual shells analysis. This is contrary to the inverse trend naively 
expected given that the Wolf Rayet stars are located at the centre of IC 132 (see On the other 
hand, it might simply mean that the WR phenomenon didn't have enough time to pollute the 
ionized gas or that, even if some contamination occurs, within the errors it cannot be observed. 
The radial increase may be explained by a) SF was propagated from the border to the central part 
of the region b) Enrichment of an old population (by shocks or SN) 'escaped' the central zone: 
and is still propagating outside. However it may be noted that (N/O) is smoothed in the cumulative 
case and also becomes compatible with a constant ratio. 

The empirical abundances in figure [26]all hint at radial increase, but again when the errors are 
considered all are compatible with a constant value through the region. All the indicators, except 
the one based on S3O3 may be set to a value of 8.0, in agreement with the mean value from the 2D 
maps. S3O3 reaches higher values by almost 0.25 dex even within the errors. 

Figure [27] shows the various line ratios used. RS2 is in almost all cases above the theoretical 
limit for the [Sn] ratio, although in the cumulative case the dispersion is lower than in the individ- 
ual one and the values are closer to the upper theoretical limit, a systematic effect that produces 

© 2012 RAS, MNRAS 000,[2^?? 



IFS ofHll regions in M 33 51 
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: Te([S III]) so. 
1 .5 - 




A B C D E F G AA AB AC AD AE AF AG 
Figure 22. IC 132: temperatures estimated for individual and accumulated shells as labeled in each panel. Vertical lines are the associated error. 

Individual Cumulative 



8 - 




12 + log(0/H) 



A B C D E F G AA AB AC AD AE AF AG 
Figure 23. IC 132: ionic and total oxygen abundances for individual and accumulated shells. 
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Individua 



Cumulative 




1.33 - 



C D E F G AA AB AC AD AE AF AG 
Figure 24. IC 132: ionic and total oxygen abundances for individual and accumulated shells. 
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Cumulative 
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A B C D E F G AA AB AC AD AE AF AG 
Figure 25. IC 132: ionic ratios and N + abundances for individual and accumulated shells. 



"unrealistic" results is evident. When comparing with the 2D maps also the trend is different as 
the spaxels at the borders of the map have low RS2 (~ 0.7). This reflects the problems with the 
line deblending, its dependence with the S/N or limitations in the atomic values to estimate n e in 
the low density regime, as discussed previously in §3.11 

Log(S0 2 3) is constant through the region; taking into account the errors the range lies within 
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A B C D E F G AA AB AC AD AE AF AG 
Figure 26. IC 132 empirical abundance indicators for individual and accumulated shells. 



the fit produced by 



Diaz et al. 



Dfaz et al. 



(2007) for T e determination. It is worth noticing that 
(|2007|) only used high metallicity Hll regions in their fit. This sets a limit log(SO 2 3)=-1.0 for 
the validity of the relation although many objects exist below this limit even in the compilation 
they present. If the fit is extended to low metallicity objects that may shift the relation to produce 
higher temperatures. The values of log(R 2 3), log(Ar 3 3 ) and log(S 3 3 ) are compatible, within the 
errors, with no gradient. The ionization sensitive ratio log(0 32 ) (McGaugh 1991) indicates that the 
ionization decreases with radius. 



by 



Figur el28lsho ws the softness parameter ?7 / =([Oii]3727/[Oiii]4959,5007)/([sii]6717,673i/[Siii]9069,9532) , defined 



Vflchez et al. 



(1988) as an indicator of the effective temperature. To first order insensitive to 
the ionization parameter in the nebula, it grows inversely proportional to the temperature, rf is also 
compatible with a constant value when errors are considered, and an increase with radius is not 
expected considering T c trends in figure |22] 

Ionization sensitive ratios are plotted in figure |29l a decrease of ionization with radius is con- 
cluded from the three indicators. 

In general when comparing the shell analysis with the 2D maps statistical agreement is found 
between the mean of the distribution and the standard deviation (s). However some evident fea- 
tures in the maps disappear in the shell analysis. A clear example is the "wall" that runs diagonally 
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Figure 27. IC 132 various line ratios for individual and accumulated shells. 
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Figure 28. IC 132 ionization parameter from sulphur and oxygen and the parameter r/' (see text) for individual and accumulated shells. 
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Figure 29. IC 132 ionization sensitive line ratios for individual and accumulated shells. 
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Figure 30. IC 132 extinction parameter C(H/3) and equivalent width EW(H/3) for individual and accumulated shells. 

in the south-east section of many 2D maps for IC 132. Hence the shell analysis must always be 
complemented with the maps for a full analysis of the region and exploitation of the IFS capabili- 
ties. 




■ T C(H/S) 


~i | £— — f — § — | - 








EW(H/S) 



7.2 Integrated properties of the Individual central regions 



In contrast to the IC 132 frame that shows a single region, our observed frames in the centre of 



M33 contain several localized Hn regions. We have used the charts from iHodge et al. 



(2002) to 



matc h separate region s in the Ha map. The identification of the regions together with the number 



from Boulesteix et al. 



(II 9741) (which we called BCLMP) is shown in Fig. I 
A segmentation analysis was performed for BCLMP 93; the rest of the central regions are 
too faint to obtain spectra with S/N high enough if segmented. The Ha levels used to measure 
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BCLMP 93 integrated shells are shown in Fig. [32] The measured line fluxes and derived quantities 
are given in tables |8l 191 [101 and ITTl while the corresponding plots are in figures |33l [341 1331 and |36l 

Figure [33] shows a steep increase of T c ([Sm]) for the individual shells but a moderate one for 
the cumulative analysis. The sulphur ionic abundances show a slight decrease with radius although 
they are consistent with being uniform when errors are considered. 

Empirical abundance indicators are represented in figure |34] N2, albeit with large errors, is 
compatible with constant abundance. The R23 based abundance decreases for the individual shells 
but becomes compatible with constant in the cumulative case. The same can be concluded for 
S3O3. Within the errors, they are all compatible with being constant. 

For the line ratios in figure [35] RS2 is always below the low density limit indicating higher 
densities in the central region than in the outskirts. They show that n e values decrease with radius 
in the individual case while for the cumulative shells the last point shows an increase in n c for the 
outermost spaxels. Even when the rest of the ratios indicate some variation trend, in general they 
are compatible with a constant value, except for the last two ratios. The log(U)[on]/[om] increases 
while r) f decreases. The decrease in if is compatible with the trend displayed by T e ([Sm]) in figure 
M 

All the ionization sensitive ratios in figure [36] show a clear increase with radius. This is puz- 
zling, as one would expect that when [Om]/H/3 increases, both [Nll]/Ha and [Sll]/Ha should 
decrease. Possible causes may be a projection effect of two zones physically separated but on the 
same line of sight or in general, multiplicity of the region or even the meddling effect of diffuse 
radiation. 

The large error bars in some line ratio plots are a consequence of the relatively small S/N of the 
data. We can nevertheless distinguish some differences between IC 132 and BCLMP 93. Firstly the 
higher Balmer decrement C(H/3) and smaller EW(H/3) in BCLMP 93 indicate higher extinction 
and either an older Hn region or, more likely, a larger contamination due to the underlying stellar 
population unrelated to the ionizing cluster in the central regions of M33, i.e. the bulge population. 
Regarding the electron temperatures, only Te([Sm])63i2 is computed for both regions and provides 
a result that suggests an outwards increasing temperature for each region. 
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Figure 31. Identification of regions from the Hodge et al. 12002) catalog. Labels are the BCLMP numbers (see text). 
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Figure 32. BCLMP 93 shell segmentation. 



Core Ha 

© 2012 RAS, MNRAS 000,|2]-?? 



Table 8. Reddening corrected fluxes for BCLMP 93 with integration over individual Ha shells. Lines are normalized to H/3 = 100 



Individual 



A 


Ion 


A 


B 


C 


D 


E 


F 


G 


H 


3727 


[O II] 


160.7±11.4 


183.9±12.6 


230.8±21.4 


232.2±18.9 


235.0±24.4 


244.8±21.9 


264.8±25.1 


294.7±52.4 


4010 


H(5 


23.4±3.2 


24.3±3.2 


21.4±4.5 


21.2±5.5 


21.0±5.1 


21.0±5.9 


17.4±7.1 


13.8±11.7 


4340 


H 7 


43.0±2.3 


43.3±2.0 


42.3±3.0 


43.0±2.9 


41.6±3.5 


43.0±4.4 


41.7±5.1 


41.3±9.5 


4363 


[0 III] 1 


-4.6±1.3 


-4.0±1.4 


-5.8±2.0 


-4.9±1.9 


-5.9±2.1 


-7.8±2.5 


-12.6±2.9 


-21.3±4.6 


4861 


[H/3] 


100.0±3.5 


100.0±3.4 


100.0±5.2 


100.0±4.7 


100.0±4.8 


100.0±5.5 


100.0±7.0 


100.0±10.8 


4959 


[0 III] 


2.4±1.0 


2.0±1.2 


3.6±1.6 


4.4±1.5 


4.9±1.3 


6.6±1.9 


9.5±2.9 


15.2±3.8 


5007 


[O III] 


9.4±0.8 


7.9±0.9 


11.1±1.2 


11.4±1.3 


13.5±1.4 


18.0±2.0 


22.1±2.6 


35.3±4.0 


6312 


[S III] 


1.3±0.9 


1.3±1.1 


1.3±1.3 


1.8±1.7 


2.1±2.1 


2.8±2.7 


3.5±3.7 


5.0±5.7 


6548 


[Nil] 


10.1±1.2 


11.0±1.5 


11.8±1.5 


13.3±1.5 


12.0±1.4 


12.2±1.2 


13.2±2.0 


12.6±3.5 


6564 


Ha 


283.8±7.6 


283.8±7.5 


281.2±10.8 


280.1±10.1 


282.3±10.1 


282.2±11.2 


284.8±14.8 


282.3±22.6 


6584 


[Nil] 


55.5±2.1 


59.0±2.4 


66.0±3.2 


68.6±3.3 


68.0±3.1 


66.8±3.1 


69.0±4.4 


68.9±6.9 


6717 


[S II] 


25.0±1.3 


30.2±1.4 


35.5±2.0 


38.2±2.2 


40.0±2.4 


41.4±2.6 


44.5±3.2 


52.1±5.3 


6731 


[S II] 


23.3±1.1 


24.3±1.2 


28.2±1.6 


29.2±1.8 


31.0±2.0 


32.4±2.1 


34.2±2.6 


38.6±4.1 


7136 


[Ar III] 


0.6±7.1 


0.7±4.6 


1.4±6.5 


1.3±9.3 


1.7±10.9 


2.2±11.9 


2.2±10.3 


8.3±13.5 


9068 


[Sin] 


27.1±14.4 


29.0±20.7 


31.7±25.3 


32.6±30.6 


32.9±32.9 


30.3±32.6 


26.8±35.1 


15.7±44.6 



I(H/3)[ergs- 1 cm- 2 10" 13 ] 21.28±0.05 18.95±0.05 18.61±0.07 18.66±0.06 18.48±0.06 18.08±0.07 18.25±0.09 18.63±0.14 

C(H/3) 0.66±0.06 0.52±0.06 0.68±0.08 0.84±0.07 1.02±0.07 1.24±0.08 1.25±0.09 1.08±0.11 

EW(H/3)[A] 38.08±6.17 46.96±6.85 47.47±6.89 55.80±7.47 50.28±7.09 48.44±6.96 42.13±6.49 33.43±5.78 



1 Contaminated by mercury street lamps. 



Table 9. Integrated temperatures, abundances and line ratios for BCLMP 93 individual Ha shells. 



Individual 





A 


B 


C 


D 


E 


F 


G 


H 


J e (ISlIl]) 63 i2 


8000±200 


0/1/1/1 1 
8000±200 


/600±300 


84UU±4UU 


n 1 / 1/ 1 1 c / 1/ 1 

91UU±5UU 


1 /n/1/1 1 h 
10300±/00 


12400±1 100 


16000±2/00 




6000±300 


6100±300 


6200±300 


6200±300 


6200±300 


6300±300 


6500±400 


6700 ±400 


12+log(S+/H+) 


6.34±0.66 


6.38±0.91 


6.50±1.05 


6.39±1.26 


6.32±1.38 


6.20±1.55 


6.05±1.95 


5.90±3.04 


12+log(S++/H+) 


6.93±0.42 


6.96±0.59 


7.04±0.68 


6.95±0.80 


6.86±0.86 


6.74±0.94 


6.54±1.11 


6.33±1.57 


12+log(0/H) (N2) 


8.45±0.8l 


O A^l 1 f\ / \ / \ 

8.47±0.90 


r i 1 j"7j j"7j ^ 

8.51±2.25 


O F **7 I A « f\ 

8.53±2.10 


( J £*/7j I -"N /~\^7 

8.52±2.07 


8.51±2.48 


8.52±4.50 


8.53±10.60 


12+l0g(0/H) (R23«pper) 


8.92±0.10 


8.87±0.l0 


8.75±0.10 


8.75±0.10 


8.75±0.10 


8.73±0.10 


8.69±0.10 


8.63±0.10 


12+log(0/H) (Ar 3 3 ) 


8.55±2.16 


8.63±0.89 


t 1 / \ ^^7 

8.67±0.57 


8.67±0.85 


O SI} 1 /-\ ^7 A 

8.68±0.74 


8.67±0.64 


C\ S' A 1 f\ S A 

8.64±0.64 


^7 r 1 r\ "i ^7 

8.75±0.13 


12+log(0/H) (S3O3) 


8.88±0.07 


8.92±0.11 


8.88±0.12 


8.87±0.l4 


O O A \ f~\ 1 A 

8.84±0.14 


AA 1 A 1 ^7 

8.80±0.13 


^7 1 t~\ 1 a 

8.75±0.l4 


8.69±0.19 


RS2 


1.07±0.01 


1.24±0.01 


1.26±0.01 


l.3l±0.0l 


1.29±0.01 


1.28±0.01 


l.30±0.0l 


1.35±0.02 


n e 


328 — 


137 — 


127 — 


85 — 


101 — 


109 — 


91 — 


55 — 


log(S0 23 ) 


-0.01±0.03 


-0.03±0.03 


-0.09±0.03 


-0.08±0.03 


-0.09±0.03 


-0.12±0.04 


-0.18±0.04 


-0.24±0.06 


N2 


-0.7l±0.ll 


-0.68±0.12 


-0.63±0.27 


-0.6l±0.25 


-0.62±0.25 


-0.63±0.30 


-0.62±0.53 


-0.61±1.25 


log(R 23 ) 


0.24±0.08 


0.29±0.07 


0.39±0.10 


0.39±0.09 


0.40±0.ll 


0.43±0.10 


f\ A ^7 1 /"\ 1 1 

0.47±0.ll 


0.54±0.l9 


log(0 3 2) 


1 1 **7 I /~\ 1 ^7 

-I.l3±0.l3 


-l.27±0.l7 


-1.20±0.16 


-I.l7±0.l5 


-l.ll±0.l5 


-l.00±0.l4 


-0.92±0.16 


/ \ ^7^7 | f\ /*7 1 

-0.77±0.2l 


log(Ar 3 3 ) 


-l.l9±3.97 


-l.03±2.ll 


-0.91 ±1.64 


-0.93±2.38 


-0.90±2.l9 


-0.9l±l.84 


-1.00±1.60 


-0.63±0.56 


log(S 3 3 ) 


0.56±0.l8 


0.65±0.26 


0.54±0.30 


0.54±0.35 


0.45±0.37 


0.32±0.4l 


0.18±0.48 


-0.05±0.68 


log(u) ([Sii]/[Siii]) 


-2.33±0.3l 


-2.39±0.43 


-2.44±0.50 


-2.47±0.59 


-2.53±0.63 


-2.57±0.68 


-2.70±0.81 


-2.85±l.l5 


log(u) ([On]/[Olll]) 


-3.93±0.05 


-4.03±0.06 


-3.98±0.06 


-3.95±0.05 


-3.9l±0.05 


-3.82±0.05 


-3.76±0.05 


-3.63±0.07 


log(V) 


l.52±0.ll 


l.63±0.l3 


1.52±0.15 


1.48±0.l6 


l.38±0.l7 


l.25±0.l7 


1.09±0.17 


0.85±0.24 


log([Oni]A5007/H/3) 


-l.03±0.04 


-l.l0±0.05 


-0.95±0.05 


-0.94±0.05 


-0.87±0.05 


-0.74±0.05 


-0.65±0.05 


-0.45 ±0.05 


log([Nn]A6584/Ha) 


-0.7l±0.0l 


-0.68±0.02 


-0.63±0.01 


-0.6l±0.02 


-0.62±0.0l 


-0.63±0.0l 


-0.62±0.02 


-0.6l±0.03 


log([Sn]AA6717,31/Ha) 


-0.77±0.07 


-0.72±0.08 


-0.64±0.08 


-0.62±0.08 


-0.60±0.08 


-0.58±0.09 


-0.56±0.09 


-0.49±0.09 



on 



Table 10. Reddening corrected fluxes for BCLMP 93 with integration over cumulative Ha shells. Lines are normalized to H/3 =100 











Cumulative 










Jesus . 






AA 


AB 


AC 


AD 


AE 


AF 


AG 


AH 


3727 


[Oil] 


160.7±11.4 


171.6±11.4 


190.3±13.8 


200.4±14.1 


207.1 ±14.7 


213.1±15.2 


220.2±16.1 


229.4±17.7 




4010 


H5 


23.4±3.2 


23.8±3.1 


23.0±3.1 


22.6±3.5 


22.3±3.7 


22.0±3.9 


21.4±4.2 


20.4±4.9 




4340 


H7 


43.0±2.3 


43.1±2.1 


42.9±2.1 


42.9±2.2 


42.6±2.3 


42.7±2.5 


42.6±2.6 


42.4±3.2 




4363 


[0 III] 1 


-4.6±1.3 


-4.3±1.3 


-4.8±1.4 


-4.8±1.4 


-5.1±1.5 


-5.5±1.6 


-6.4±1.7 


-8.1±2.0 


re 
3 


4861 


[H/3] 


100.0±3.5 


100.0±3.4 


100.0±3.7 


100.0±3.9 


100.0±3.9 


100.0±4.1 


100.0±4.4 


100.0±5.0 


4959 


[O III] 


2.4±1.0 


2.2±1.1 


2.7±1.1 


3.1±1.2 


3.4±1.2 


3.9±1.2 


4.7±1.4 


6.0±1.6 




5007 


[O III] 


9.4±0.8 


8.7±0.8 


q s^n q 

y. JILlYJ.y 


Q Q-t-H Q 
y.yztzu.y 


10.6±1.0 


U.8±l.l 


13.2±1.2 


15.9±1.3 


idez 


6312 


[S III] 




1 . jm 1 .u 


1 3^ 1 1 


1 A4-1 9 




1 7+1 (\ 


z.uzn L.y 


9 4.-1-9 A 


6548 


[Nil] 


in 1 -t-1 9 


1U. JZtZ 1 - J 


11.1±1.3 


11.4±1.3 


I 1 'X 

I I .JZtZ 1 .J 


1 1 <C_1_1 9 


1 1 Q_l_ 1 
1 l.Oll.J 


1 1 Q_l_ 1 < 


re 


6564 


Ha 


283.8±7.6 


283.8±7.4 


283.0±8.0 


282.3±8.3 


282.3±8.4 


282.3±8.6 


282.6±9.2 


282.6±10.5 


J 4 


6584 


[Nil] 


55.5±2.1 


57.1 ±2.2 


60.0±2.4 


62.1±2.5 


63.2±2.5 


63.8±2.5 


64.5±2.7 


65.0±3.1 




6717 


[S II] 


25.7±1.3 


28.8±1.3 


30.6±1.4 


32.3±1.5 


33.6±1.6 


34.7±1.7 


34.6±1.8 


38.2±2.1 




6731 


tsm 


23.7±1.1 


23.1±1.1 


24.2±1.1 


25.7±1.2 


26.7±1.3 


27.0±1.4 


29.9±1.5 


29.5±1.7 




7136 


[Ar III] 


0.6±7.1 


0.7±5.6 


1.0±3.8 


1.0±3.8 


1.1 ±4.7 


1.3±5.6 


1.4±6.1 


1.7±6.9 




9068 


[Sm] 


27.1±14.4 


27.7±17.2 


28.8±19.5 


29.9±22.0 


30.7±24.0 


30.7±25.1 


30.1±26.5 


29.8±29.1 




ICH^Cergs-^m- 2 10" 13 ] 


21.28±0.05 


40.23±0.10 


58.83±0.16 


77.49±0.21 


95.97±0.27 


114.05±0.33 


132.30±0.41 


150.92±0.53 




C(H/3) 




0.66±0.06 


0.60±0.04 


0.62±0.04 


0.68±0.03 


0.75±0.03 


0.83±0.03 


0.89±0.03 


0.92±0.03 




EW(H/3)[A] 


38.08±6.17 


41.80±6.47 


43.44±6.59 


45.89±6.77 


46.67±6.83 


46.94±6.85 


46.21±6.80 


44.13±6.64 





1 Contaminated by mercury street lamps. 




* -0 

* -0 



Table 11. Integrated temperatures, abundances and line ratios for BCLMP 93 accumulated Ha shells. 



Cumulative 





AA 


AB 


AC 


AD 


AE 


AF 


AG 


AH 


xursni/n^io 

J eVL VJAA1 J/Dolz 


8000±200 

Uuvv L ^— V./ V./ 


8000±200 

Uvuv L » V./ V./ 


7900+200 


8000+300 


8200+300 


8600+300 


9100+400 

y 1 v./ v / 1 r v./ v./ 


9900+500 


Te([SlIl])s0 23 


6000±300 


6000±300 


6100+300 


6100+300 


6200+300 


6200+300 


6200+300 


6300+300 


12+log(S+/H+) 


6.35±0.66 


6.35±0.77 


6.39+0.86 


6.39+0.96 


6.37+1.03 


6.33+1.10 


6.29+1.21 


6.21 + 1.36 


12+log(S++/H+) 


6.93±0.42 


6.95±0.50 


6.97+0.56 


6.96+0.62 


6.94+0.66 


6.90+0.70 


6.84+0.75 


6.76+0.83 


J. z. i ivj ti V vjv li-/ > z.y 


O .T\J 1 W. < ] J. 


o 4A+0 82 

O .TAJ J_U.Oi. 


o 40+1 05 


8.49+1.19 


8.49+1.25 


c 50+1 11 


c 50+1 56 

O --IM 1 J. .^J VJ 


c 50+2 06 

O . -J VJ _1_Z- . VJVJ 


LZ.^ivg\\Ji ri) yr^23upper ) 


c Q2+0 1 
o . ~z. i w. 1 vj 


c qn+0 1 


c 85+0 1 

O. O.J 1 w. 1 V/ 


81+0 1 

O.OJXu. 1 VJ 


CI +0 1 

O.Ol 1 w. 1 VJ 


80+0 10 

O . uu_Lu. 1U 


7Q+0 1 

0. / y 1 v / . 1 vj 


77+O 1 

O. / / _1_VJ. 1 VJ 




o 55+? 16 

O . _1_Z. . J. V.J 


8 61+1 18 

O . VJ J. I I.IO 


f.4+0 51 


f.4+0 51 


A5+0 55 

O . VJ^y I V /- / 


fi5+o 56 

. vj.j 1 >>.-») 


A5+0 57 

O . VJ.J _1_VJ.»7 / 


A5+0 52 

O . VJ,J _1_VJ . Z- 




8 88+0 07 


o Qfl+n 09 

o . y v.* i v./ . v/ y 


8Q+0 10 
. y _l_vj. j. vj 


O 8Q+0 1 1 
. oy 1 * ' - J. J. 


8 88+0 1 1 

O.OO _I_VJ. J. J. 


8A+0 12 


8,4+0 12 
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Figure 33. BCLMP 93 temperature and ionic abundances estimated from the sulphur lines for the individual and cumulative shells. 

Due to the radial decrease in the ionization parameter, we expect to observe a radial decrease in 
ionization; IC 132 (fig. [29]) shows a decreasing [Olll]/H/3 and increasing [Nll]/Ha and [SII]/Ha, 
but in BCLMP 93 (fig.l36l) these line ratios seem to increase outwards. We have already discussed 
a possible explanation due to confusion in the central region. 

8 CONCLUSIONS 

The main objective of this work was to study with spatially resolved spectroscopy the differences 
in the characteristics between the innermost and outermost regions of star formation in nearby 
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Figure 34. BCLMP 93 empirical abundance estimates for the individual and cumulative shells. 
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Figure 35. BCLMP 93 line ratios for the individual and cumulative shells. 

spiral galaxies with the aim of exploring the differences in the properties of intense star formation 
activity in widely contrasting environments. 

We have presented the results of our IFS data obtained with PMAS in PPak mode at the CAHA 
3.5m telescope, of Hn regions in M33, selected to cover a wide range in environment and proper- 
ties of the star forming regions. Our sample includes the luminous Hn region IC 132 in the outer 
disc and the innermost Hn regions in a field of 500pc x 300pc that surrounds the galaxy centre. 
The results are presented in the form of intensity maps of the Hn regions intrinsic parameters: line 
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Figure 36. BCLMP 93 ionization ratios for the individual and cumulative shells. 
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Figure 37. BCLMP 93 extinction parameter and EW(H/3) for the individual and cumulative shells. 

ratios, excitation, extinction, electron density, electron temperatures and elemental abundances 
with a sampling of l"x I". We also studied internal gradients of the intrinsic parameters of the Hn 
regions for apertures or shells defined according to their average surface brightness in Ha. 
Our main findings are: 

Extinction 

We find a clear difference in the extinction and its radial gradient between the innermost and outer 
regions in M 33. While IC 132 shows a global value of C(H/3) of about 0.3 and indications of lower 
extinction at larger radius, C(H/3) for BCLMP 93 is globally around 1.3, decreasing outwards in 
the region. 
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In both systems the EW(H/3) is smaller in the centre of the region and grows to a maximum 
in the intermediate shells before dropping in the outer shells. The radial behaviour of the EW(H/?) 
can be understood by a combination of factors, firstly the continuum in the inner regions is strongly 
affected by the massive ionizing stars, secondly the outer regions are affected by the disc (IC 132) 
or bulge (BCLMP 93) component and globally the brightness of the emission lines decreases 
radially outwards. Thus the EW(H/3) is probably peaking in the transition zone between the edge 
of the ionizing cluster and the regions where the surface brightness of H/3 drops below that of the 
older population. 

We find a smaller equivalent width of H/3 in BCLMP 93 than in IC 132 indicating higher 
extinction and either an older Hll region or more likely a larger contamination due to the presence 
of an old stellar population in the central regions of M 33, i.e. the bulge population. 

It is interesting to note that the maximum EW(H /?) ~ 300 A in IC 132 is close to the maximum 



of ~ 400 A for SB99 models (Leitherer et al. 



1999), assuming instantaneous star formation with 



upper mass limit = 1OOM , zero age main sequence, an initial mass function with a = 2.35 and 
Z=0.020Z Q . The asymptotic maximum of the cumulative EW(H/3) reaches ~ 200 A indicating an 
ionizing cluster with an age of less than 3 Myrs that is in agreement with the detection of WN stars 
that are expected to appear after 2.5 Myrs of evolution of a young massive cluster. 

For BCLMP 93 the same procedure suggests an upper limit for the age of the cluster of about 
5 Myrs, bearing in mind that the large contamination due to the underlying bulge population pre- 
cludes a more precise age estimate. This age makes the central region older than IC 132 and could 
explain the lack of present day WR features in the region. 



Electron density 

The electron density (n e ) was estimated from the RS2=[Sll]A6717/[Sll]A6731 line ratio after de- 
blending the doublet. This diagnostic situates the Hll regions in the low density limit. Even if 
an exact value for n c is not possible, a variation in the ratio with position is observed, support- 
ing the fact that to assume a single n c value to characterize a whole Hll region is too simplistic. 
Better spectral resolution, access to alternative density diagnostics and revision of the theoretical 
estimation based on new atomic data for collisional strengths are necessary to shed light on the 
subject. 
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Electron temperature 

For the central zone of IC 132 the electron temperature (T e ) was estimated from the detected 
[Om]A4363 and [Sm]A6312 lines allowing the use of the [Om]AA4959,5007/[Om]A4363 and 
[Sill] AA9069,9532/[Sm] A6013 ratios to compute T e ([OlIl]) and Te([Sm]) respectively. For BCLMP 93 
the lack of detection of [Olll] A4363 implies that only Te([Slll]) was determined in its central zone. 
The blueshift of the galaxy means that the [0m] A 4363 line is only an upper limit, as it is affected 
by the street light mercury line (Hg I A 4358.3). 

The 2D distribution of T e ([Olll]) and T e ([Oll]) in IC 132 appears to be rather uniform in the 
high surface brightness region surrounded by slightly hotter gas. The mean values are 16100 K for 
T e ([Olll]) and between 15300 K and 13700 K for T e ([Oll]) with an rms scatter per spaxel of about 
3200 K in all cases. For T e ([Slll]) we find a mean value of 16200 K with an rms of 3800 K over 
279 spaxels. Given that the maps include about 300 spaxels, the error of the mean temperatures 
is less than 600 K. These results are confirmed by the segmentation analysis that shows a positive 
gradient in all temperature estimates and a global cumulative value of 13500 K for T e ([0n]) and 
15100 for T e ([Olll]) while T c ([Sm]) reaches an asymptotic value of 13600 K. 

Higher T c ([Sm]) is derived for IC 132 than for BCLMP 93. This is in agreement with the 
centre being more metallic than the disc. Consequently in central regions cooling is more efficient 
through fine structure lines in the infrared. 

The analysis in shells indicates that the direct method is reliable for the brightest zones. How- 
ever as the radial distance within the regions increases the lines are still detected in the integrated 
spectra and can be included in the usual temperature computation but the inclusion of fainter zones 
produces physically unrealistic temperature values. 

This can be seen in tables [7] and [TT] and in figures [22] and [33] where the outer ring drives the 
temperatures up. This is true for temperature determinations using [Olll]A4363 and [Oll]A6312 
and hence can not be attributed just to the contamination of [Olll]A4363 by street lights due to the 
particular blueshift of M33. 

Abundances 

The smallest dispersion in O/H determined from empirical indicators is obtained for 03N2 fol- 
lowed by N2. Although the related line ratio may have a large dispersion along the region, whether 
or not it is transferred to (O/H) depends on the calibration of the relation. If the conversion is lin- 
eal, the steepness of the slope dictates the dispersion in (O/H). This is why the oxygen abundances 
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from 03N2 and N2 are almost constant despite the fact that the l ine ratios have bee n discarded as 



2010). 



good indicators because their larger dependence on excitation ARelano et al.\ 

R23 displays a dispersion of ~ 0.3 dex for the central zone, however when individual spax- 
els are converted to (O/H) an extra difficulty is found because the upper and lower branches are 
included for the same object, thus the spaxels in the tra nsition zone hav e larger uncertainty de- 
pending on the calibration chosen. This is worse for the IPilyuginl (12001 ah method because a first 
estimation of (O/H) is required and the area of transition between branches is not defined. The fact 
that for the same region both branches are involved is suspicious and raises the question about the 
validity of the method when applied to individual spaxels. 

Sulphur is the only direct indicator in common for IC 132 and BCLMP 93. It shows higher 
abundance for the central region coinciding with the results for the rest of the empirical indicators. 



Diagnostic diagrams 

Our search for systematic differences between the inner and outer Hn regions in M33 was com- 
plemented with the use of diagnostic diagrams. In particular we focused in the [Olll]A 5007/H/?, 
vs. [Nil] A 6584/Ha diagram and found that the individual spaxels display a wide range of ion- 
ization conditions, extending by more than 2.5 dex from the location of the brightest spot. The 
distribution of points for the two Hn regions considered displays a striking clear dichotomy. While 
the points for BCLMP 93 show a vertical scatter, the points for IC 132 have a horizontal spread 
orthogonal to that of BCLMP 93. The regions occupied by the distributions have hardly any super- 
position while some points clearly move away from the Hn region zone of the diagnostic diagram 
intruding in the AGN and transition zone. This intrusion raises a question about the wisdom of 
using the results of global models in the diagnostic diagrams for individual spaxels of resolved 
systems, including the definition of the HII boundary line. 

Thus it is important to be cautious when applying or comparing model predictions for global 
observations to spatially resolved data. In this case the BPT diagnostic diagrams may give mis- 
eading results for sp axels with ionization conditions that may not be reached in global spectra 



( Ercolano et al. 



2012). 



Bearing this in mind we computed pseudo 3D photoionization models and compared the dis- 
tributions of observed line ratios in the diagnostic diagram with the results of our photoionization 
models and found that it was possible to reproduce the spread observed in the [Olll]5007/H/3 ver- 
sus [Nll]6584/Ha diagram for both zones. We conclude that the position and horizontal spread of 
points in IC 132 is reproduced by models with high U, high T* and low Z while the position and 
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the vertical scatter of the observed points in the central Hll regions is better reproduced by models 
with lower U and and higher Z. The contribution of the central X-ray source to the ionization 
budget is not clear, since not very anomalous line ratios were detected in its proximity. 

Wolf Rayet stars 

Blue WR features were isolated and characterized for IC 132, the spatial resolution is enough 
to distinguish two concentrations. No WR features were detected in the central zone in apparent 
discordance with theoretical predictions that higher metallicity regions should have larger WR 
content, but perhaps consistent with the estimated larger age of the central region. 
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